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ABSTRACT 

We present a systematic analysis of 43 nearby galaxy groups (fcTsoo = 0.7 — 2.7 keV or A/500 = 
10^^ - 10^''/i~^ M0, 0.012 < z < 0.12), based on Chandra archival data. With robust background 
subtraction and modeling, we trace gas properties to at least r25oo for all 43 groups. For 11 groups, 
gas properties can be robustly derived to rsoo- For an additional 12 groups, we derive gas properties 
to at least riooo and estimate properties at rsoo from extrapolation. We show that in spite of the large 
variation in temperature profiles inside 0.15 rsoo, the temperature profiles of these groups are similar at 
> 0.15 rsoo and are consistent with a "universal temperature profile." We present the K — T relations 
at six characteristic radii (30 kpc, 0.15 rsoo, ^2500, ''isooj ''1000 and r5oo), for 43 groups from this work 
and 14 clusters from the Vikhlinin et al. (2008) sample. Despite large scatter in the entropy values at 30 
kpc and 0.15 rsoo, the intrinsic scatter at r25oo is much smaller and remains the same (~ 10%) to rsoo- 
The entropy excess at rsoo is confirmed, in both groups and clusters, but the magnitude is smaller than 
previous ROSAT and ASCA results. We also present scaling relations for the gas fraction. It appears 
that the average gas fraction between r25oo and rsoo has no temperature dependence, ~ 0.12 for 1-10 
keV systems. The group gas fractions within r25oo are generally low and have large scatter. This work 
shows that the difference of groups from hotter clusters stems from the difficulty of compressing group 
gas inside of r25oo ■ The large scatter of the group gas fraction within r25oo causes large scatter in the 
group entropy around the center and may be responsible for the large scatter of the group luminosities. 
Nevertheless, the groups appear more regular and more like clusters beyond r25oo, from the results on gas 
fraction and entropy. Therefore, mass proxies can be extended into low mass systems. The M500 — Isoo 
and M500 — ix 500 relations derived in this work are indeed well behaved down to at least 2 xlO^"^ h~^ 

Subject headings: cosmology: observations — dark matter — galaxies: clusters: general — X-rays: 
galaxies: clusters 



1. INTRODUCTION 

Most baryons in clusters reside in the hot ICM, and most 
clusters are low-mass groups and poor clusters because 
the mass function of virialized systems is steep. Studies 
of galaxy groups are thus especially important for under- 
standing the gravitational and thermal evolution of the 
bulk of matter in the Universe. In contrast to hot clus- 
ters, galaxy groups are systems where baryon physics (e.g., 
cooling, galactic winds and AGN feedback) begins to dom- 
inate over gravity. Groups are not just scaled-down version 
of massive clusters. Cluster scaling relations (e.g., L — T 
and entropy - T) show deviations from the self-similar rela- 
tions at the low-mass end (reviewed in Voit 2005), refiect- 
ing the importance of baryon physics, which is essential 
to shape the properties of the cluster galaxies and ICM. 
Cooling sets a characteristic entropy threshold in the ICM 
(Voit & Ponman 2003) and is required to understand the 
cluster scaling relations (Voit 2005). On the other hand, 
cooling has to be balanced by feedback to preserve the ob- 
served cluster stellar fraction and the galaxy luminosity 
function. The most important feedback may be the AGN 
outflows driven by the central SMBH. AGN outflows can 
simultaneously explain the antihierarchical quenching of 
star formation in massive galaxies, the exponential cut-off 
at the bright end of the galaxy luminosity function, and 
the quenching of cooling-flows in cluster cores (e.g., Scan- 
napieco et al. 2005; Croton et al. 2006). Along with 



SN winds, they act to suppress cooling and structure for- 
mation. They pump thermal energy into the surrounding 
ICM as the host galaxy and the central SMBH formed and 
evolve. These imprints are recorded in the ICM and are 
reflected in the cluster scaling relations. 

There are two important ICM properties that directly 
reflect the role of baryon physics. The flrst is ICM en- 
tropy With ROSAT and ASCA data, Ponman et al. 
(2003) showed that at 0.1 r2oo, the ICM entropy (K) 
obeys a simple relation with the cluster temperature (T), 
K (X T~0-65^ which is different from the self-similar rela- 
tion {K oc T). Ponman et al. (2003) and Voit & Ponman 
(2003) demonstrated that this K — T relation at 0.1 r2oo 
can be understood with a simple model involving a fixed 
cooling threshold, which may be related to galaxy forma- 
tion. However, the details of cooling and feedback that 
determine this threshold and its scatter are still unclear. 
Finoguenov et al. (2002) and Ponman et al. (2003) showed 
that groups and clusters have significant excess entropy at 
'"500- Voit & Ponman (2003) suggested that smooth ac- 
cretion (via e.g., galactic winds in subhalos) can boost en- 
tropy higher than lumpy hierarchical accretion to produce 
the observed excess. However, Borgani et al. (2005) show 
that galactic winds alone are not able to boost entropy 
very much at radii beyond rsoo. The AGN-like feedback 
algorithm in Kay (2004) has a more substantial effect on 
entropy at large radii, which may imply that radio outflows 
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have an important impact on group properties. Thus, it is 
important to robustly constrain the entropy around r5oo 
with observations and to connect the dispersion with other 
group properties. The Ponman et al. (2003) results rely 
on extrapolations based on assumed density (/3-model) and 
temperature (polytropic-model) profiles, which may not be 
adequate and may bias the derived ICM properties (e.g., 
Borgani et al. 2004; Vikhlinin et al. 2006, V06 hereafter). 
These early results can be examined with Chandra and 
y>MM.- Newton. There has been some work on the ICM en- 
tropy of groups with the Chandra or XMM-A''ewio77. data 
(e.g., Mushotzky et al. 2003; Sun et al. 2003; Pratt et al. 
2005, 2006; Mahdavi et al. 2005; Finoguenov et al. 2007). 
Our much larger sample allows us to study the scatter in 
group properties. As the ICM properties are traced to at 
least r25oo for each group, the K — T relation at <; r25oo 
can be better constrained from our studies. 

Another important ICM property is the gas fraction. 
Groups have low gas fr-actions within r25oo (0.03-0.07, V06; 
Gastaldello et al. 2007, G07 hereafter). However, few 
groups have gas fraction measured to r5oo. V06 and G07 
derived gas fraction within rsoo for four low-temperature 
systems {kT < 2.7 keV). The results, 0.06 - 0.15, span a 
large range. Gas fraction in simulated groups and clusters 
is directly related to the strength of cooling and star forma- 
tion (e.g., Kravtsov et al. 2005), and a small gas fraction 
in groups may imply efficient cooling and star formation. 
The enclosed gas fraction can also be modified by the AGN 
feedback (e.g., Puchwein et al. 2008), thus bearing the im- 
print of the feedback history. It is therefore of great inter- 
est to know how gas fractions and total baryon fractions 
(gas -I- stars) of groups at rsoo compare with those of clus- 
ters. The fgas — T relation is also important for cosmology, 
e.g., determining cosmological parameters, and using the 
Sunyaev-Zel'dovich flux as a proxy for the cluster mass 
(e.g., Vikhlinin et al. 2008, V08 hereafter). Besides the 
science related to entropy and gas fraction, better deter- 
mination of the group gas properties are also essential to 
constraining the low-mass ends of other important scaling 
relations (e.g., M - T,M - Y-^ and c - M). M - T and 
M — Ix relations are essential for using clusters to study 
cosmology (e.g., Kravtsov et al. 2006; V08), while the 
c — M relation is important for understanding the forma- 
tion and evolution of the dark matter halos (e.g., Buote et 
al. 2007). 

We started a systematic analysis of galaxy groups in the 
Chandra archive to better constrain the ICM properties 
in groups and to better understand the difference between 
groups and clusters. In this paper, we present the results 
on 43 galaxy groups that appear relaxed beyond the cen- 
tral core. Our sample is larger than in the previous work 
on 3 - 16 groups with Chandra and XMM-A^ewton data 
(Mahdavi et al. 2005; Pratt et al. 2005; Finoguenov et 
al. 2007; G07; Rasmussen & Ponman 2007). The data 
are homogeneously analyzed with results on temperature, 
entropy, gas mass and total mass. The sample size also al- 
lows us to measure the scatter in various interesting ICM 
properties. We focus on the data reduction and the derived 
scaling relations in this paper, while more extensive dis- 
cussions with modeling and work on an extended sample 
including irregular groups will be presented in future pa- 
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pers. The group sample is defined in §2. The data analysis 
is presented in §3, including spatial and spectral analysis. 
We especially discuss the Chandra background and our 
method of background subtraction in this section and Ap- 
pendix. In §4, we define four tiers of groups with different 
data coverage. Different characteristic radii (r2500j '^isooj 
riooo and rsoo) are reached in different tiers. We also de- 
fine the system temperatures in §4 (Tsoo etc.) and derive 
their empirical relations. The group temperature profiles 
are discussed in §5, while §6 is about the ICM entropy. In 
§7, we discuss M — T and M — Ix relations, gas fraction, 
concentration parameter, baryon fraction and fossil groups 
in this sample. There are groups with signs of AGN heat- 
ing and groups with strong central radio sources in this 
sample. We discuss them in §8. Systematic errors are 
discussed in §9. §10 is the summary and conclusions. We 
assumed Hq = 73 km s~i Mpc^\ f^M=0.24, and Oa=0.76. 

2. THE GROUP SAMPLE 

The groups in the sample are mainly selected from the 
Chandra archive as of September, 2008. We also include 
proprietary data on three groups from the Chandra pro- 
gram 09800349 (PI: Vikhlinin). We started to collect 
groups from several previous group samples constructed 
from ROSAT observations: Mahdavi et al. (2000), 
Mulchaey et al. (2003) and GEMS (Osmond, & Ponman 
2004). We have also searched for low temperature sys- 
tems (T < 3 keV) in the on-line database BAX0. However, 
most ROSAT groups with well-determined temperatures 
are quite nearby. The most distant group in GEMS is at 
z=0.0282 (NGC 6338). As we want to constrain gas prop- 
erties out to at least r25oo (ideally rsoo) for each system, 
many nearby systems in these ROSAT catalogs are ex- 
cluded. Thus, we have had to extensively examine the data 
in the Chandra archive seeking groups with good observa- 
tions out to these large radii. Our selection criteria are: 1) 
full-frame Advanced CCD Imaging Spectrometer (ACIS) 
data at the focal plane temperature of -120C (after Jan. 
29, 2000); 2) T500 ^ 2.7 keV and a global temperature con- 
strained to better than 15%; 3) 0.015 < z < 0.13 and group 
emission traced to at least r25oo with the Chandra data; 4) 
group emission well centered around the cD and not sig- 
nificantly elongated or disturbed beyond the group core. 
T500 is the temperature measured between 0.15 r5oo and 
''500 (defined and discussed in §4) . The upper limit on Tsoo 
is determined from the M500 — T500 relation in V06 to make 
sure that we are studying low-mass systems with M500 of 
< lQ^'^h~'^ M0. We understand that there is not a well- 
defined temperature boundary separating galaxy groups 
and clusters. Many people may consider groups to be sys- 
tems with temperatures of < 2 keV. Nevertheless, we refer 
to all the systems in our sample as groups for convenient 
comparison with the clusters in V08. The requirement to 
well constrain T500 allows us to derive a temperature pro- 
file for each group. The constraint on redshift enables a 
first cut to make sure r25oo can be reached. We used the 
r — T relation in V06 as the first guess. It turns out the de- 
rived r~T relations from our work are close to V06's. With 
the r25oo — T relation in V06, r25oo — 11.6' for an 1 keV 
group at z = 0.015. An ACIS pointing can reach this ra- 
dius although the coverage is partial. The only exception 
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to the redshift requirement is NGC 1550 at z = 0.0124. 
There are two offset observations for this luminous group 
that allow us to reach riooo- Clearly some hotter systems 
at z ^ 0.015 may not meet our criterion to reach r25oo and 
are excluded after their temperatures are constrained. We 
indeed have examined many more systems than the ones 
in our final sample. With Chandra^s superior angular res- 
olution, almost all groups have substructures around the 
center, and in some cases these features can be very sig- 
nificant and striking, like the central 60 kpc of IC 1262 
(~ 0.1r5oo). Nevertheless, they are included as long as the 
group emission is well centered and regular at ^ 0.15r5oo, 
where our main interests are. Our science goals concern- 
ing gas entropy are also not affected by the dynamic state 
of the group. Moreover, at least the Fx mass proxy is still 
robust in unrelaxed systems (Kravtsov et al. 2006). 

The Chandra archival search not only includes targets 
in the cluster and galaxy categories, but also targets in 
the AGN category, as many radio galaxies are in group 
environments. There are also targets selected optically. 
Many faint systems are included, as shown by the wide 
range of gas entropy at 30 kpc and 0.15 r^oo (§6). Thus, 
our final sample is not much biased to the X-ray luminous 
systems. The final sample includes 43 groups, listed in 
Table 1. It is hardly a homogeneous or a complete sam- 
ple. But the wide spread of the gas densities and entropy 
values at r < 0.15r5oo indicates a wide luminosity range 
as cool cores in groups contribute a large portion of their 
luminosities (See §7.2). If there were groups that are much 
fainter or poorer in gas than the faint systems in this sam- 
ple, it is difficult to constrain their gas properties with the 
current X-ray instruments anyway. Because of the layout 
of ACIS, full coverage at large radii (defined as > 80% cov- 
erage for the outermost spectral extraction bin, as there is 
always area missing from extended and point sources and 
chip gaps) is only achieved in 7 groups (Table 2). However, 
the coverage at r ^ ^2500 is much better. Previous analy- 
sis for nearby clusters or groups with the Chandra data all 
had partial coverage at large radii (e.g., V06; Rasmussen 
et al. 2007; V08) and our sample is about 3 times bigger 
than these previous ones. 

3. THE DATA ANALYSIS 

3.1. Chandra observations & calibration steps 

All observations were performed with ACIS. Standard 
data analysis was performed which includes the corrections 
for the slow gain change^ and Charge Transfer Inefficiency 
(for both the FI and BI chips). We investigated the light 
curve of source-free regions (or regions with a small frac- 
tion of the source emission) to identify and exclude time 
intervals with particle background flares, including weak 
flares. The relevant information on the Chandra point- 
ings is listed in Table 1. We corrected for the ACIS low 
energy quantum efficiency (QE) degradation due to the 
contamination on ACIS's Optical Blocking Filter 0, which 
increases with time and is positionally dependent. The 
dead area effect on the FI chips, caused by cosmic rays, 
has also been corrected. As the background subtraction is 
important for this project, we present it in detail in the 

* http:/ /cxc. harvard. odu/contrib/alexey/tgain/tgain. html 
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next section and Appendix. We do not use any data on 
the S4 chip, because of the residual streaks often seen at 
low levels (after running the CIAO tool "DESTREAK") 
and the lack of the stowed background data. The "readout 
artifact" , seen in groups with very bright cores, is also cor- 
rected (see e.g., V06). We used CIA03.4 for the data anal- 
ysis. The calibration files used correspond to Chandra Cal- 
ibration Database (CALDB) 3.4.3 from the Chandra X-ray 
Center, released in March, 2008. We are aware of a pos- 
sible over-correction of the Chandra effective area beyond 
the Iridium M-edge (~ 2 keV) in CALDB 3.4.3 and before 
(see David's presentation in the 2007 Chandra Calibra- 
tion Workshop) 0, which can bias Chandra temperatures 
to higher values, especially for hot clusters (T > 4 keV). 
The difference between Chandra and XMM-Newton tem- 
peratures is also shown in Snowden et al. (2008). The cali- 
bration work to incorporate this into the Chandra CALDB 
is ongoing. However, as shown in both David's presenta- 
tion and Snowden et al. (2008), the agreement between 
Chandra and XMM-Newton temperatures for ^ 4-5 keV 
systems is very good. The highest gas temperature at any 
radii in our sample is ~ 3.3 keV for several groups at the 
center, while temperatures at large radii are much lower 
and mainly determined by the iron L hump (instead of the 
continuum slope as for hot clusters). Thus, any changes 
in our results from this correction should be smaller than 
the current measurement errors and scatter. In the spec- 
tral analysis, a lower energy cut of 0.4 keV (for the BI 
data) and 0.5 keV (for the FI data) is used to minimize 
the effects of calibration uncertainties at low energy. The 
solar photospheric abundance table by Anders & Grevesse 
(1989) is used in the spectral fits. Uncertainties quoted in 
this paper are 1 a. 

3.2. Determining the local background 

Proper background subtraction is important for the 
analysis in the low surface brightness regions of the groups. 
A detailed discussion of the Chandra background and the 
relevant dataset is presented in the Appendix. We de- 
termine the local background based on the stowed back- 
ground data and modeling. The corresponding stowed 
background of each observation is scaled according to the 
flux ratio in the 9.5 - 12 keV band (e.g., V05). We adopt a 
3% uncertainty on the particle background (PB) normal- 
ization (5% for BI data in period E, see Appendix) in the 
error budgets. 

The local cosmic X-ray background (CXB) can be mod- 
eled. The best fit of the hard CXB component is deter- 
mined with an absorbed power law with a photon index 
of 1.5. The absorption is determined from spectral fits to 
the group spectra, which involves iterations and is present 
in §3.3. The soft CXB component is adequately described 
by two thermal components at zero redshift, one unab- 
sorbed component with a fixed temperature of 0.1 keV 
and another absorbed component with a temperature ei- 
ther derived from spectral fits or fixed at 0.25 keV (see 
Appendix). Abundances of both components are fixed at 
Solar. We can compare our CXB model with other work. 
V05 used the blank sky background and corrected for dif- 
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ferences between a local control field and the blank sky 
background. As the scaling factors of the PB are close to 
unity in their sample, V05 ignored the correction for the 
unresolved hard CXB. However, the correction for hard 
CXB becomes more important after the middle of 2004 
when the PB flux is 30% - 50% higher than that in the 
blank sky background (Appendix). V05 found that gener- 
ally a single thermal component with a temperature of ~ 
0.2 keV and solar abundance can fit the soft CXB excess 
or decrement well. In regions with high RASS R4-5 values, 
a second thermal component with a temperature of ~ 0.4 
keV and solar abundance is required. Humphrey & Buote 
(2006) and G07 also used two thermal components with 
solar abundances to describe the soft CXB. The CXB tem- 
peratures are fixed at 0.07 kcV and 0.2 kcV. As discussed 
in the Appendix, Snowden et al. (2008) used three thermal 
components to describe the soft CXB. Thus, our model of 
the soft CXB only has less freedom than the model by 
Snowden et al. (2008). However, our model balances the 
requirement of having enough components to fit the local 
CXB and avoiding parameter degeneracy in generally low- 
statistics data. First, the two ~ 0.1 keV components will 
naturally be mixed in our energy band when the absorp- 
tion is low. Second, statistics for most of our groups at 
large radii are not very good, so uncertainties are mainly 
statistical errors. As long as the background uncertainties 
are folded into the final error budgets, very detailed mod- 
cling of the soft CXB component is not crucial. Third, 
we test different models of the soft CXB for groups in 
our sample for which a large group-emission-free region is 
available. All spectra can be well fit by our two-component 
model, while the three-component model by Snowden et 
al. (2008) makes little or no improvement. In fact, the soft 
CXB can be well fitted by a single thermal component in 
most cases. However, the abundance usually has to be free 
and the best-fit value is usually very close to zero. There- 
fore, we conclude that the two-component soft CXB model 
is adequate and also necessary for our analysis. 

The Chandra observations in our sample can be classi- 
fied into two categories, one with regions free of group 
emission in the FOV, another with group emission de- 
tected to the edge of the field. Cool systems at higher 
redshifts (e.g., z > 0.04) usually have group-emission- free 
regions in the off-center chips (e.g., S2 for ACIS-I obser- 
vations, 12/3 and SI for ACIS-S observations). There are 
also two groups with a(ljac;ent Chandra pointings for un- 
related targets in coincidence (ESQ 306-017 and A1692). 
We took the following approach to look for group-emission- 
free regions. The radial surface brightness profile is derived 
with the exposure correction and background subtraction 
using the scaled stowed background. Different pointings 
for the same group are combined. The obtained surface 
brightness profile has group emission plus the CXB. The 
region where CXB is more dominant than the group emis- 
sion, if present, can be determined from the flattened por- 
tion at the outer region of the surface brightness profile. 
We used a power law plus a constant at large radii to de- 
termine whether a significant group-emission-free region 
exists and the radial range of that region. If such a region 
is found, we have a control field with only CXB. In this 
work, the inner radius of the group-emission-free region is 
> 3' larger than the outermost radius for the group tem- 



perature profile. This control field may be a single region 
on the S2 chip (for the ACIS-I data) or two separate re- 
gions on the SI and 12/13 chips respectively (for the ACIS- 
S data). Although we could simply fit the spectrum (or 
spectra) of this control field to determine the local CXB, 
the statistics are generally not sufficient. Moreover, the 
hard CXB in this control field may be larger than that 
in the outermost bin for the group temperature profile, as 
the flux of the hard CXB depends on point source exci- 
sion. Thus, we fit the spectra of the control field and the 
outermost bin for the group temperature profile together 
to better constrain the local CXB. The normalizations of 
the soft CXB in different regions are linked by the ratios 
of their covered solid angles of the sky, while the normal- 
izations of the hard CXB are not linked. Therefore, we 
can determine the local CXB in the outermost bin for the 
group temperature profile. In this work, we assume that 
the soft CXB is constant across the examined group area. 
Accurate determination of the local CXB is only impor- 
tant for large radii. The covered area at large radii tends 
to be in the same direction from the group center so the 
assumption should be reasonable. There is a complication 
that the hard CXB may be smaller in inner radial bins but 
it is a small effect and can be corrected (also discussed later 
in this section). 

However, many groups in our sample have group emis- 
sion detected to the edge of the field, so we have to fit the 
CXB components with the group emission together. Since 
these groups are limited by the Chandra FOV, the group 
emission is generally still significant near the edge of the 
field. We fit the spectra from the two outermost radial bins 
together, with the normalizations of the soft CXB compo- 
nents linked. Generally in this work, we are conservative 
and do not include the outermost bin in the temperature 
profile as the uncertainties are generally large. However, 
there are a few groups that remain bright to the edge of 
the field (e.g., NGC 1550, A262 and MKW4). They all 
have SI data, so the soft CXB component can be easily 
separated from the group emission (e.g., Fig. 1), owing 
to the good response of the BI chips at < 1 keV. As the 
temperatures of the cool gas can be well constrained from 
the iron hump centroid, group emission can be robustly 
separated from the soft CXB even to the edge of the field. 
Thus, we derive temperature profiles to the edge of the 
field for these groups. 

As wc! separate the PB from the CXB, the absolute CXB 
flux can be derived in each group field. In Table 2, we list 
the flux of the local background components, for both the 
soft CXB and the hard unresolved CXB (in the outermost 
bin). We expect the derived soft CXB is on average cor- 
related with the RASS R45 flux in the surrounding area 
(excluding the source region, see Table 2), which is true as 
shown in Fig. 2. One should be aware of the uncertainty 
in the R45 flux, as the RASS soft X-ray background maps 
have 0.2 deg pixels and the various uncertainties combined 
are not small (e.g., variable SWCX emission). There are 
also uncertainties related to cross-calibration and cosmic 
variance of both the soft and hard CXB, so a detailed one- 
to-one comparison is hardly meaningful. The derived 2- 8 
keV flux of the hard CXB component depends on the lim- 
iting flux of the observation. As discussed in Appendix, 
there is an empirical relation between the limiting flux for 
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point sources and tlie average 2-8 IceV flux of unresolved 
hard CXB (K07). We estimated the average limiting flux 
of the group observations in the outermost spectral bin. 
The regions we used to control the local background gen- 
erally have an area of about one ACIS chip. As shown 
in Hickox & Markevitch 2006 (HM06 hereafter) and K07, 
the cosmic variance in this angular scale (depending on the 
two-point angular correlation function of point sources) is 
20% - 30% and the total hard CXB flux has ~ 10% un- 
certainty. Thus, the expected hard CXB flux from the 
empirical relation is only meaningful in an average sense. 
We use the CIAO tool MKPSF to generate several PSFs 
(at an energy of 1.4 keV) in the outermost spectral bin. We 
then derive the 90% enclosed power aperture and measure 
the 3(7 limits at these regions. Their average is taken as the 
estimate for the limiting flux for point sources. Compared 
with the average growth curve determined by K07 (their 
Fig. 19), we have an estimate of the unresolved hard CXB, 
which is also listed in Table 2. We can sec from Table 2 
that the general agreement is quite good, while the uncer- 
tainties for the hard X-ray CXB are much larger than the 
difference. It is also true that the limiting flux for point 
sources depends on the off-axis angle. The Chandra lim- 
iting flux changes little within the central 6-7' from the 
aimpoint, but increases rapidly beyond 10'. This can be 
corrected from the slope of limiting flux vs. unresolved 
hard CXB in K07 and the absolute flux in the outermost 
bin. The correction is small, also because the errors of the 
hard CXB are not small. 

Besides the 3% uncertainty on the PB mentioned above 
(5% for BI data in period E), we also included the following 
error budgets and added them (and the PB uncertainty) to 
the statistical uncertainties of temperatures, abundances 
and surface brightness in quadrature. For the hard CXB, 
we set the photon index at 1.1 and 1.9 (e.g., HM06) and 
repeat the analysis to estimate the local background in 
those situations. For the soft CXB component, if the tem- 
perature of the hotter component cannot be derived from 
the data and is fixed at 0.25 keV, we change it to 0.2 keV 
and 0.3 keV and repeat the analysis to estimate the local 
background in those situations. Thus, we are conservative 
in the uncertainties for the local CXB. 

3.3. Spectral analysis and deprojection 

Once the local CXB is determined, wc proceed to de- 
rive the projected temperature profile. The radial bins for 
temperature measurement are decided from the outermost 
bin, by requiring the temperature to be constrained to bet- 
ter than 30% (1 a with all error budgets) and Tout I fin = 
1.3 - 1.75. After determining the outermost bin, the inner 
bins are determined progressively with Tout I Tin = 1-25 - 
1.6. Point sources and sub-clumps are excluded. The ab- 
sorption column is determined from the MEKAL fit to the 
integrated spectrum between O.lrsoo and 0.4r5oo. In this 
process, we also tried the VMEKAL model with extra free 
parameters of O, Ne, Mg, Si, S, Fe and Ni abundances, as 
a lower 0/Fe ratio may cause excess absorption (Buote et 
al. 2003). This process and the determination of the local 
CXB are done in iterations to obtain the final value of the 
absorption. If the derived absorption is consistent with the 
Galactic absorption from the HI survey (Table 2) within 
1 C7, the Galactic value is used. Only when both MEKAL 



and VMEKAL fits show excess absorption, the X-ray ab- 
sorption from the MEKAL fit (always consistent with the 
VMEKAL fit within 1 a) is adopted (Table 2). About 44% 
of groups show excess absorption relative to the Galac- 
tic value and indeed we find IRAS 100 fiia enhancement 
in many of these groups (see Appendix for discussions of 
some groups). We also examined the absorption variation 
with radius in each group. Beyond the central 20 kpc, 
no significant absorption variation is found for any group 
so we used a fixed absorption column for each group. In 
several cases, we observed an absorption increase within 
the central 20 kpc, but this analysis is complicated by the 
possible multi-phase gas around the center. Nevertheless, 
it has little effect on any of our results. We also discuss 
the systematic error related to absorption in §9. 

We used the MEKAL model to fit the spectra of the 
group emission. The free parameters are temperature, 
abundance and normalization, once the absorption is de- 
termined. For spectra around the center, we also include 
a component to account for the LMXB emission from the 
cD galaxy. The LMXB component is represented by a 
power law with an index of 1.7, within D25 aperture ob- 
tained from HypcrLeda. The total LMXB luminosity is 
fixed from the Lx — -^Ks relation derived in Kim & Fab- 
biano (2004), where Lks is the total Kg band luminosity 
of the cD from 2MASS. We also assume that the LMXB 
emission follows the Kg band light (see also Gilfanov 2004). 

The projected group temperature profiles are shown in 
Fig. 3-5. We used the algorithm determined by Vikhlinin 
(2006) to derive the deprojected temperature profiles in 
a parametric way, which was first applied in V06. The 
required inputs are the three-dimensional (3D) or depro- 
jected profiles of gas density and abundance, and the pro- 
jected temperature profile. V06 simply used the projected 
abundance profile as their sample is dominated by hot 
clusters. However, the emissivity of ^ 2 keV plasma 
is sensitive to the chemical abundances so the 3D abun- 
dance profile is required for our work. We applied the 
non-parametric geometrical deprojection (summarized in 
e.g., Pizzolato et al. 2003; G07) to derive deprojected 
abundances in wider radial bins, generally merging 2-3 
adjacent bins for the temperature profile. The deprojected 
abundance is an emission-weighted average in each bin so 
an effective radius in each bin is required. We define the 
effective radius as the emissivity-weighted radius in each 
bin, where the plasma emissivity, e{T,Z,r), depends on 
the 3D profiles of temperature, abundance and density. 

Since the determination of the 3D profiles of tempera- 
ture, abundance and density depends on each other, iter- 
ations have to be done to derive the best-fit 3D profiles. 
Usually at most three iterations are required to stabilize 
the best-fits of these profiles. We assumed the following 
3D abundance profile: 

Z{r) = Zo + Ziexpi-ir/rz)"^) (!) 

This simple form can fit the abundance profiles of all 
groups in the sample. An example is shown in Fig. 6 
that represents the best-constrained 3D abundance profile 
in this sample (thus the most difiicult to fit, as the errors 
are the smallest). Most groups in this sample only have 
3D abundances constrained in 3 - 5 bins with larger er- 
rors so good fits can always be achieved with this simple 
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function. Once the best-fits of all profiles are achieved, we 
applied 1000 Monte Carlo simulations to address the un- 
certainties of the 3D abundance profile. The simulation is 
realized by scattering the 3D abundance profile according 
to the measurement errors. In this process, to be conser- 
vative, we also include an uncertainty of 10% of the bin 
size on the effective radius of each bin for the abundance 
profile. Thus, besides the best-fit 3D abundance profile, 
we have 1000 simulated profiles to cover the error ranges, 
which will be used to estimate the uncertainties of the 3D 
temperature and density profiles. 

We used the same form of the 3D temperature profile as 
used in V06: 

T(r\ - ^o(^Acooi)°-' + Trnin (r/n)-° 

The exceptions are A1139, A1238 and A2092, which are 
faintest 2-3 keV systems in this sample so the errors on 
temperatures are large. Their 3D temperature profiles are 
modeled as: T(r) = To[l + (r/rt)'']^/", which fits their 
temperature profiles very well (Fig. 4 and 5). The inner 
boundary of the fit is 5 - 20 kpc, depending on the qual- 
ity of the fits at the center. As we focus on properties at 
large radii, the detailed choice of the inner boundary (5 
- 20 kpc) is not a concern. Because the sizes of bins at 
large radii are not small, each bin is further divided into 
6 sub bins. The temperature modeling is done at the cen- 
ter of each sub bin, only added up (with the algorithm by 
Vikhlinin 2006) at the end to obtain the expected tem- 
perature for each original bin of spectral analysis. In this 
way, the binsize effect is minimized and we don't need to 
worry about the accurate determination of effective radius 
for each radial bin, as most previous work had to address. 
The uncertainty of the 3D temperature profile is derived 
from 1000 Monte Carlo simulations. In each simulation, 
the measured projected temperature profile is scattered ac- 
cording to the measurement errors, and a new simulated 
abundance profile is input. The density profile is fixed at 
the best-fit value, as the temperature error is the dominant 
error source. The reconstructed 3D temperature profiles, 
with 1 a uncertainties from simulations, are also shown in 
Fig. 3-5. There are groups where a central corona exists 
(e.g.. Sun et al. 2007), so naturally the temperature gradi- 
ent at 5 - 20 kpc is large in these cases (e.g., A2462, A160, 
ESO 306-107, HOG 51 and NGC 6269). Sometimes the 
3D temperature appears too low within 5-10 kpc radius 
(e.g., A160 and NGC 6269), which however affects none of 
our results in this work as core properties are excluded. 

3.4. Gas density 

We extract the surface brightness profile in the 0.7 - 
2 keV band, as suggested by V06 to avoid the 0.6 keV 
hump in the soft CXB. Point sources and chip gaps are 
excluded. The scaled stowed background is subtracted so 
only the group emission and the local CXB is left. Some 
previous analysis on a surface brightness profile involved 
the correction of a single exposure map. However, the 
Chandra exposure map is energy and position dependent. 
It is not accurate to use a single exposure map (even one 
convolved with the group spectrum) as there is spectral 
variation in a group. The most accurate approach is to 



generate response files for each bin of the surface bright- 
ness and convert the raw count rate (without corrections 
on vignetting and other response) to density, from the 
derived 3D temperature and abundance profiles and re- 
sponse files. To achieve that, we use XSPEC to generate 
an MEKAL emissivity matrix that depends on tempera- 
ture, abundance and position of the radial bin for surface 
brightness (or response files there), once the absorption is 
determined. The ranges for temperatures and abundances 
cover the observed ranges for any particular group. This 
emissivity matrix provides the conversion factor needed to 
transfer the observed count rate to the emission measure. 
We assume the following density model (with 11 free pa- 
rameters) : 

„2 _ ^2 (r/rc)-" 1 

° [1 + (r/r,)2]3/3-a/2 [1 + {r/r,)^]^h 

+«02[l-(V^c2)"1'^^/"=(r <re2) (3) 
This model combines the profile used in V06 and the 
profile proposed by Ettori (2000) for cool cores. We find 
that this model provides very good fits for all groups in our 
sample. The density profile is then converted to the emis- 
sivity profile from the emissivity matrix. The emissivity 
profile is then projected along line of sight with the for- 
mula of geometrical deprojection (e.g., McLaughlin 1999). 
The local CXB can be added later with the known CXB 
spectra and the radial set of the response files. The result- 
ing surface brightness profile can be compared with the 
observed one. To avoid the binsize effect that may espe- 
cially affect outskirts as wider bins are required there, we 
further divide bins within 100 kpc to 3 sub bins and bins 
beyond the central 100 kpc to 6 sub bins. The conversion 
and deprojection is done in these sub bins. They are later 
merged to compare with the observed profile. The density 
errors are estimated from 1000 Monte Carlo simulations, 
with 1000 corresponding simulated abundance and tem- 
perature profiles as inputs. 

This method is similar to that used in V06, but we 
use the 3D temperature and abundance profiles to con- 
vert count rates to density. Like V06, we also analyzed 
the ROSAT PSPC pointed observations for the purpose 
of constraining gas density profiles at large radii. The in- 
clusion of the PSPC data is especially useful to z < 0.04 
groups. We used the software developed by Snowden et 
al. (1994) to produce fiat-fielded PSPC images in the 0.7 
- 2 keV band (more exactly, R567 bands, or PI channels of 
70 - 201). The images are further analyzed as described in 
Vikhlinin et al. (1999). Seventeen groups have sufficiently 
long PSPC data (listed in Table 1 with effective exposure) 
and we included the PSPC surface brightness profile in the 
modeling of the density profile. The probed outermost ra- 
dius of the PSPC data for each group is listed in Table 2. 
We only used the PSPC surface brightness profiles outside 
the central 3' to avoid the PSF correction in the core. In 
all cases, there is good agreement between the Chandra 
and PSPC surface brightness data. One example is shown 
in Fig. 7. 

4. CHARACTERISTIC RADII AND DEFINITION OF SYSTEM 
TEMPERATURES 

As stated above, we generated 1000 simulated profiles 
of the 3D temperature, abundance and density profiles. 
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which cover the ranges of the measurement errors. As- 
suming hydrostatic equihbrium, each set of the 3D tem- 
perature and density profiles determine a set of character- 
istic radii of the group {r^ao, ''looo, ''1500 and r25oo)- The 
profiles of interesting quantities, e.g., total mass, gas mass 
and entropy, and their values at the characteristic radii 
are also determined. For each quantity, the peak in the 
distribution from 1000 simulations defines the most prob- 
able value (or "best fit"). The la errors at two sides are 
estimated by determining the regions that contain 68% of 
realizations at each side. The results are listed in Table 3 
and 4. r5oo is the basic characteristic radius that is used 
in most scaling relations in this work. However, it cannot 
be robustly determined for all groups in our sample so the 
^"500 — T scaling relation is required. As the radial coverage 
of data in this sample differs from one group to another, we 
define four tiers of groups for which group properties (en- 
tropy, mass, Yx, and gas fraction) are derived to different 
characteristic radius. 

• Tier 1: groups with surface brightness (including 
the PSPC profile) derived at > 2a levels to > r5oo 
(note we are conservative to estimate the errors of 
the local CXB, see notes of Table 2) and a tem- 
perature profile derived to > 80% of r^oo- Eleven 
groups in our sample are in this tier and their tem- 
perature profiles are derived to 81% - 117% of r5oo 
with a median of 97%. Groups in this tier have den- 
sity profiles derived to > t^qq so the density gradi- 
ent at ~ r5oo is constrained. Temperature profiles 
are derived beyond r^oo or sufficiently close to rsoo 
for reasonable extrapolation. We note that in V06, 
there are three clusters and groups (A383, MKW4 
and A 1991) with temperature profiles derived only 
to 73% - 89% of rsoo, but properties at rsoo are still 
derived from extrapolation. Thus, our criteria are 
similar to those in V06 and V08. 

• Tier 2: groups with surface brightness and temper- 
ature profiles all derived to at least ripoo but not in 
tier 1. Twelve groups are in this tier and their tem- 
perature profiles are derived to 68% - 87% of rsoo 
with a median of 77%. We consider that groups in 
this tier have properties well constrained close to 
rsoo so rsoo is determined in these groups. Group 
properties at rsoo are also derived from extrapola- 
tion. However, in figures of scaling relations at rsoo, 
tier 2 groups are always marked differently from tier 
1 groups. The fits with or without them are both 
listed. 

• Tier 3: groups with surface brightness and temper- 
ature profiles derived to at least risoo but not in 
tiers 1 and 2. Eleven groups are in this tier. Nine 
of them have temperature profiles derived to 52% 
- 72% of rsoo, which is close to riooo (~ 0.73rsoo 
in this sample). Thus, group properties at riooo 
arc also derived from extrapolation. A1238 and 
RXJ 1206-0744 are put in this tier as the temper- 
ature and density errors are large although riooo is 
reached in both cases. 

• Tier 4: nine other groups with surface brightness 
and temperature profiles derived to at least r2soo- 



Group properties at risoo are also derived from ex- 
trapolation in this tier. 

Previous X-ray work on clusters and groups often de- 
fined < r >, which is the emission- weighted temperature 
within a certain aperture. As definitions are generally dif- 
ferent, it is necessary to use an unified definition that is 
easily accessible from observations. In this work, we define 
the system temperature as: 

• Tsoo: the spectroscopic temperature measured from 
the integrated spectrum in the projected 0.15 rsoo 
- rsoo annulus. 

We derive Tsoo by integrating the 3D temperature pro- 
file from 0.15 rsoo to 1.6 rsoo (or riso), in an annular 
cylinder with projected radii of 0.15 rsoo - '''soo along the 
line of sight, with the algorithm by Vikhlinin (2006). The 
choice of the outer radius in the 3D integration little affects 
Tsoo as it is emission- weighted. The same definition of the 
system temperature was also used in Nagai, Kravtsov & 
Vikhlinin (2007b, NKV07 hereafter), Maughan (2007) and 
V08. This definition excludes the central region, where a 
cool core or a locally heated region may exist. Indeed 
the group temperature profiles are much more similar at 
r > 0.15rsoo (§5 and Fig. 8). This temperature can also be 
directly derived from data, provided that the full coverage 
of rsoo is achieved (generally not the case for groups in our 
sample). Similarly, we can define Tiooo, ^isoo and T2soo, 
with the projected inner boundary always at 0.15 rsoo and 
the projected outer boundary at riooo, ''1500 and r25oo re- 
spectively (3D temperature profile still integrated to 1.6 
rsoo)- Their empirical relations can also be determined. 

As we cannot derive rsoo for groups in tiers 3 and 4, the 
'"500 — Tsoo relation needs to be determined. The r — T 
relation is just a manifestation of the M — T relation so it 
is presented in §7.1. We also derived the average ratios of 
the characteristic radii from 23 tier 1 and 2 groups: riooo / 
rsoo = 0.741±0.013, nsoo / rsoo = 0.617±0.011 and rasoo 
/ rsoo = 0.471±0.009, which are about what are expected 
for the average C500 of this sample (^ 4.2, §7.3), 0.727, 
0.599 and 0.465 respectively (assuming an NFW profile). 

We also need empirical relations between Tsoo and Tisoo 
(or T2500) to estimate Tsoo for groups in tiers 3 and 4. For 
23 groups in the first and second tiers, we found: 

fcTsoo/fcTisoo = 0.93 ± 0.02 (4) 

This ratio is not temperature dependent in our sample 
and the fit is very good (x^/dof = 6.6/22). If we only 
fit 11 groups in the first tier, the ratio is the same. We 
notice that V08 derived a similar empirical relation be- 
tween Tsoo and temperature measured at 0.15 - 0.5 rsoo 
(close to our T25oo)- V08 also included a linear term as 
they explored a wider temperature range. We also derived 
fcT5oo/fcT25oo = 0.89 ± 0.02, but the scatter is larger as 
shown by the poorer fit (x^/dof = 18.6/22). The empiri- 
cal relation between fcTsoo and fcTisoo, combined with the 
rsoo ^ ^500 scaling, allows us to estimate rsoo for groups in 
tiers 3 and 4, in a few iterations. We specifically used the 
rsoo~2500 relation from 23 tier 1 and 2 groups and 14 clus- 
ters in V06 and V08 (The fourth row of Table 6). For refer- 
ence, we also give the best fit of the r2soo — ^500 relation for 
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all 43 groups in this sample: {E{z)r25QQ /155 ± 4h ^kpc) = 
(T5oo/lkeV)° '520±o.o4o^ ^^^^^ jg gimiiar to the VOG result, 
(i5(^)r250o/146±3h-ikpc) = (T5oo/lkeV)0-547±o.o20_ 

5. TEMPERATURE PROFILES 

Scaled temperature profiles for these groups are shown 
in Fig. 8, in logarithmic and linear scales. We scale tem- 
peratures with T2500, which is robustly determined for each 
group. While the scatter within the central 0.15 rsoo is 
large, the group temperature profiles are more similar be- 
yond 0.15 rsoo, with a declining slope similar to that pre- 
dicted from simulations (e.g., Loken et al. 2002). From 
0.15 rsoo to ~ rsoo, we can fit the projected temperature 
profiles with this simple form: 

T/T2S00 = (1.22 ± 0.02) - (0.79 ± 0.04)r/r5oo (5) 

Interestingly, V05 derived an average form of T/ < 
T >= 1.22 - 1.2r/ri8o (0.125 < r/nso < 0.6) for 13 
systems in their sample, which is very similar to ours as 
rsoo ^ 0.62ri8o. If we use a similar form as used in Loken 
et al. (2002), we can also fit the projected temperature 
profiles at 0.15 rsoo - fsoo with this form: 

T/r25oo = (1.37 ± 0.03) (1 + r/rsoo)~('-''^°-''' (6) 

Based on ASCA data, Markevitch et al. (1998) first 
suggested that temperature profiles of clusters are self- 
similarly declining with radius. This result was later con- 
firmed by De Grandi & Molendi (2002) with the Bep- 
poSAX data, and by V05 with the Chandra data, and by 
Piffaretti et al. (2005), Pratt et al. (2007) and Leccardi 
& Molendi (2008, LM08 hereafter) with the XMM-Newton 
data. Self-similarly declining temperature profiles are also 
generally observed in simulations (e.g., Loken et al. 2002; 
Borgani et al. 2004; Kay et al. 2004). The self-similar 
decline in groups was also suggested by Sun et al. (2003), 
G07 and Rasmussen & Ponman (2007). Our group sample 
is the largest one so far with detailed studies. As shown 
in Fig. 8, group temperature profiles are generally self- 
similar with a slope consistent with that in simulations, 
although there is some scatter. We also combine all data 
points to make a mean temperature profile (Fig. 9). Fig. 
9 also plots the mean temperature profile of hot clusters 
at z=0. 1-0.3 from LM08, as well as the mean tempera- 
ture profile of 1-3 keV systems from the Borgani et al. 
(2004) simulations. We do not plot the mean tempera- 
ture profiles from De Grandi & Molendi (2002), V05 and 
Pratt et al. (2007), but they are close to LM08's (Fig. 
21 in LM08). Clearly, the group temperature profiles are 
more peaky than those of clusters, starting to decline at 
~ 0.2rsooj where the mean cluster temperature profile is 
still flat. This difference was also noticed by V06. The 
group temperature profiles actually agree more with sim- 
ulations (e.g., Loken et al. 2002; Borgani et al. 2004), 
which generally have problems explaining the flat temper- 
ature profiles of clusters at 0.15r5oo — 0.3r5oo (e.g., Borgani 
et al. 2004). 

6. ENTROPY 

The ICM entropy, defined as K = T/rie^^, is a fun- 
damentally important quality of the ICM as summarized 



in Voit (2005). Entropy records the thermodynamic his- 
tory of the ICM, as non-gravitational processes (e.g., AGN 
heating, cooling and star formation) deviate the entropy 
relations from the self-similar relations determined with 
only gravity and accretion shocks (e.g., Ponman et al. 
1999; Voit & Bryan 2001). Voit et al. (2005, VKB05 
hereafter) derived the baseline ICM entropy profile from 
simulations in the absence of non-gravitational processes, 
K{r)/K2oo.Mi = 1.32(r/r2oo)^-^- ^200, adi is an entropy 
scale of a non-radiative cluster at r2oo (VKB05). This 
baseline entropy profile can be compared with our obser- 
vational results to measure the impact of non-gravitational 
processes on the ICM. With a weighted average of ^ 4.2 
for csoo in this sample (§7.3), r5oo/r2oo = 0.669 for an 
NEW mass profile and the baseline entropy relation is con- 
verted to: i^(r)/i^5oo,adi = l-40(r/r5oo)^'^ (note that we 
use /Tsoo.adi here for the adiabatic entropy scale at rsoo, 
see the definition of -ft'200.adi in VKB05 and Voit 2005). 
For the cosmology assumed in this paper and a baryon 
fraction (/b) of 0.165 (Komatsu et al. 2008), 

irsoo,adi = 342 l,eYcm^{^^r/'Eiz)-'/%,'^' (7) 
Note that -ft'soo.adi has a different h dependence from the 

— 1/3 

observed entropy {K oc h^^ ). From the derived temper- 
ature and density profiles, we obtained entropy profiles for 
each group. The best-fit scaled entropy profiles are shown 
in Fig. 10. Substantial scatter is present within 0.3 - 0.4 
rsoo- As Fig. 10 is basically a plot showing the large scat- 
ter of the density profiles, it is clear that there is a wide 
range of X-ray luminosities in this sample, which is also 
implied by the K — T relations shown next. We discuss 
the ICM entropy of these groups in the following three 
sections. 

6.1. K — T relations 

We first examined the K — T relations at 
30 kpc, 0.15 rsoo, ''2500, ''isoo, ?'iooo and rsoo 
(Fig. 11). Entropy values at these radii are: 
^30kpc, -fi'o.isrsoo, ^2500, ^1500 , i^iooo , -f^'soo respectively 
(note the difference between ifsoo and -ftr5oo,adi defined 
earlier). We chose 30 kpc to represent the entropy level 
around the BCG. Previous studies often discussed entropy 
at 0.1 r2oo (e.g., Ponman et al. 2003), which is about 0.15 
'"500- We also performed the BCES (Y|X) regression (Akri- 
tas & Bershady 1996) to these relations and the results 
are listed in Table 5. The entropy values at 30 kpc radius 
hardly show any correlation with the system temperature. 
The K — T relation is stronger at 0.15 rsoo, but the in- 
trinsic scatter is still substantial. The large scatter of gas 
entropy around the center has been known for a while 
(e.g., Ponman et al. 2003). We also examined the con- 
nection between the luminosity of the central radio source 
(see Table 1) and the entropy scatter at 0.15 rsoo- No 
correlation is found, before or after the temperature de- 
pendence of entropy is removed. However in §7.2, we show 
that the entropy scatter at 0.15 rsoo is tightly correlated 
with the scatter of the gas fraction within r25oo- 

Despite large entropy scatter from the core to at least 
0.15 rsoo, the K — T relations are tighter at r25oo and be- 
yond. We derived the intrinsic scatter in these relations, 
using the method described in Pratt et al. (2006) and 
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Maughan (2007). As shown in Fig. 11, the intrinsic scat- 
ter reduces to 10% at r25oo and stays at that level to rsoo- 
As the derived intrinsic scatter decreases with increasing 
measurement errors, could this trend be due to the in- 
creasing measurement errors with radius? To answer this 
question, we compared ^^2500 - T2500 and i^o.isrsoo - ^2500 
relations. If we use the relative measurement errors of 
-f''o.i5r50o for the corresponding ii'2500, the intrinsic scatter 
only increases to 12%. Similarly, if we use the relative mea- 
sm'cmcnt errors of i^2500 for the corresponding -fi'o.iSrSOOi 
the intrinsic scatter only decreases to 28%. Therefore, the 
significant tightening of the K — T relation from 0.15 r5oo 
to r25oo is not caused by different measurement errors. 
Similarly, we find that the intrinsic scatter in entropy at 
radii from r25oo to riooo is affected little by switching mea- 
surement errors. The K^qq values have large measurement 
errors and half of systems are extrapolated from riooo- 
Clearly more groups with -ftTsoo robustly determined are re- 
quired, but the intrinsic scatter in the current i^soo — T500 
relation is consistent with the level from r25oo to riooo- 
Thus, groups behave more regularly from r25oo outward 
than inside 0.15 rsoo- 

We also include 14 clusters from V06 and V08 for com- 
parison and to constrain the K — T relations in a wider 
temperature range (Fig. 12 and Table 5). The measured 
slopes increase from 0.50 at 0.15 rsoo to ~ 1 at rsoo, 
mainly caused by the excess entropy of groups at their 
centers. The slope we find at 0.15 rsoo, 0.494±0.047, is 
consistent with what Pratt et al. (2006) found at 0.1 
r2oo for 10 groups and clusters with the XMM-Newton 
data, 0.49±0.15, but smaller than what Ponman et al. 
(2003) found at 0.1 r2oo with ROSAT and ASCA data, 
~ 0.65. We however caution that the V08 cluster sample 
lacks non-cool-core clusters. Pratt et al. (2006) also gave 
the K — T relation at 0.3 r200) which is about r25oo- The 
slope they foimd, 0.64±0.11, is also consistent with our 
result, 0.740±0.027. At rsoo, the slope we found is con- 
sistent with the expected value from the self-similar rela- 
tion (1.0). Interestingly, the derived iiTsoo — T500 relation 
agrees very well with the NKV07 simulations with cooling 
-|- star-formation (Fig. 12), although the agreement is pro- 
gressively worse with decreasing radius. The tightening of 
the K — T relation at r25oo and beyond was also reported 
in the NKV07 simulations, but the predicted /r2500 — 7" re- 
lation by NKV07 lies below all groups in our sample (Fig. 
11 and 12). The NKV07 sinndations with cooling + star- 
formation achieve entropy amplification with strong con- 
densation to drop dense materials out of the X-ray phase. 
The resulting stellar fraction is about twice the observed 
value as too much material has cooled. If the star for- 
mation is suppressed with more efficient feedback, entropy 
is lower in the NKV07 simulations, between the results 
from the simulations with cooling -|- star-formation and 
the non-radiative simulations, which further disagrees with 
observations. It seems that the most challenging task is to 
explain the excess entropy of groups at r25oo- As shown 
in Borgani & Viel (2008), simply increasing entropy floor 
from pre-heating produces too large voids in the Lyman-a 
forest. Thus, it is still an open question on how to generate 
enough entropy in groups from the center to rsoo and still 
preserve other relations like the condensed baryon fraction 
and properties of the Lyman-a forest. 



6.2. Entropy ratios 

We derived the ratios of the observed entropy to the 
baseline entropy from VKB05 at rsoO; riooo ^-^d r2soo (Fig. 
13). For comparison, we also include the entropy ratios for 
14 clusters from V06 and V08. The entropy ratios are al- 
ways larger than or comparable to unity at all radii. The 
average ratio decreases with radius, for both groups and 
clusters, although the decrease is more rapid in groups. 
For groups, the weighted mean decreases from 2.2 at r25oo 
to 1.57 - 1.60 at rsoo; while the mean decreases only ~ 
9% for clusters. Thus, there is still a significant entropy 
excess over the VKB05 entropy baseline at rsoo, in both 
groups and clusters. The weighted mean entropy of groups 
at rsoo is ~ 18% larger than that of clusters. If we con- 
sider only the Tsoo ^1.4 keV groups, the weighted mean 
for 7 groups is 1.68. Similar studies have been done with 
ASCA and ROSAT data (Finoguenov et al. 2002; Ponman 
et al. 2003). When the observed entropy values at rsoo 

are scaled with Mc^^J^ (which is proportional to -ftT^g adi)' 
both works found that groups have on average twice the 
scaled entropy of clusters. Ponman et al. (2003) further 
concluded that the excess entropy is observed in the full 
mass range. Thus, the results from this work, V06 and 
V08 confirm the previous finding of the excess entropy at 
''500 over the full mass range, but the magnitude of the 
excess is smaller and the difference between groups and 
clusters is also smaller than previous results, ^ 15% - 20% 
vs. ^ 100% in Finoguenov et al. (2002) and Ponman et 
al. (2003). 

As pointed out by Pratt et al. (2006), the magnitude 
of the excess may be aff'ected by the systematics of the 
M — T relation, or how robust the hydrostatic equilibrium 
(HSE) mass is. The NKV07 simulations show that the 
HSE mass is systematically lower than the real mass and 
the difference is biggest in low-mass systems. The real 
mass is 45% higher than the HSE mass for 1 keV groups, 
while the difference reduces to 13% for 10 keV clusters. 
The best-fit Msoo — T500 relation from this work and V06 
(§7.1) is close to the relation for the HSE mass in NKV07. 
Thus, if this bias is real, the actual entropy ratios at rsoo 
are smaller. For 1 keV groups, a 45% higher total mass 
means a 13% larger rsoo- For a typical entropy slope of 0.7 
in this sample (§6.3), the entropy ratio at rsoo is reduced 
by 17%. Similarly for 5 keV clusters, the entropy ratio at 
rsoo is reduced by 8% for an average entropy slope of 0.9 
at large radii. Thus, this bias can explain only part of the 
entropy excess observed. The entropy baseline we adopted 
in this work is from the SPH simulations, while the AMR 
simulations produce a baseline with 7% higher normal- 
ization (VKB05) , likely because of its better capability to 
catch shocks. However, the entropy excess for groups at 
rsoo remains significant. 

There are several mechanisms to achieve entropy am- 
plification at rsoo- The first idea relies on modification of 
accretion, without extra non-gravitational processes (Pon- 
man et al. 2003; Voit et al. 2003; Voit & Ponman 2003). 
If preheating or feedback in small subhalos that are being 
accreted can eject gas out of the halo and thicken the fil- 
aments significantly, the accretion may be smoother than 
the lumpy accretion in hierarchical mergers. Voit (2005) 
showed that the smooth accretion can generate ~ 50% 
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more entropy throughout the cluster than would lumpy 
hierarchical accretion. However, Borgani et al. (2005) 
argued that this entropy amplification effect is substan- 
tially reduced by cooling. The other ideas resort to non- 
gravitational processes. As discussed in the last section, 
although the NKV07 simulations predict the i^soo — T500 
relation very well, they produce too many stars. Thus, 
models with only cooling may not be enough. Borgani 
et al. (2005) showed that galactic winds from SN explo- 
sions arc rather localized and cannot boost entropy enough 
at large radii. Thus, a feedback mechanism that can dis- 
tribute heat in a very diffuse way is required. As groups are 
smaller than hot clusters, AGN outflows from the central 
galaxy can reach a larger scaled radius. One good exam- 
ple in this sample is 3C449 (§8 and Appendix), with radio 
lobes extending to at least 3.7 r^oo- Thus, it is still an open 
question to explain the entropy excess at r5oo, especially 
in groups. But the much reduced entropy ratios from this 
work largely alleviate the problem. It is also clear that 
more systems with the entropy ratio constrained at t^qq 
are required for better comparison, especially T500 1.4 
keV groups. 

6.3. Entropy slopes 

We also derived the entropy slopes at 30 kpc— 0.15r5oo, 
0.15 r5oo-r250o, '"2500-'ri500 (orrdct.spo) and ri5oo-rdot,spe 
(Fig. 14). The scatter is large but the average slopes are 
about the same (~ 0.7) beyond 0.15 r5oo- The slope is al- 
ways shallower than that from pure gravitational processes 
(~ 1.1). Mahdavi et al. (2005) analyzed the XMM- TVewion 
data of 8 nearby groups and found their entropy profiles 
are best fitted by a broken power law with the break ra- 
dius of ~ O.lrsoo- Across the break radius, the entropy 
slope decreases from 0.92 to 0.42. As the break radius is 
very near the core, this kind of entropy profiles may be 
caused by a cool core within 0.1 tsoo- From 0.15 rsoo, 
there are systems in our sample with an entropy slope of 
around 0.42, but the weighted mean is significantly larger. 
The measured average slopes in this work are consistent 
with the result by Finoguenov et al. (2007) (0.6 - 0.7) for 
groups. 

7. MASS AND GAS FRACTION 

7.1. M -T and M -Yx relations 

One of the most important aspects of cluster science is 
to use clusters to study cosmology, which often involves 
derivation of a cluster mass function. Assuming hydro- 
static equilibrium. X-ray observations can be used to de- 
rive the cluster mass (at least the HSE value), provided 
the radial coverage of the data is good. Often we esti- 
mate cluster mass from another X-ray observable used as 
a mass proxy. A frequently used mass proxy is gas temper- 
ature so it is important to understand the cluster M — T 
relation. There had been a lot of work on the cluster 
M ~T relation before the Chandra and XMM-Newton era 
(e.g., Finoguenov et al. 2001; Sanderson et al. 2003). 
However, as emphasized in V06, it is crucial to constrain 
the gas properties (both the temperature gradient and the 
density gradient) at large radii (e.g., around r^oa)- Sys- 
tematics like the assumption of a polytropic equation of 
state and the inadequate fit of the density profile at large 
radii can bias the HSE mass to lower values (Borgani et 



al. 2004; V06). The M -T relations have also been con- 
strained with Chandra and XMM-iVewtofi (Arnaud et al. 
2005 with the XMM-Newton data; V06 and V08 with the 
Chandra data). However, the number of clusters used to 
constrain the M — T relation is still small (10 in Arnaud 
et al. 2005 and 17 in V08). There are only 4 systems with 
temperatures of 2.1 - 3.0 keV in Arnaud et al. (2005) and 
only three systems with temperatures of 1.6 - 3.0 keV in 
V08 (two overlapping with the Arnaud et al. sample). Our 
paper adds more < 2.7 kcV systems (23 tier 1/2 groups). 

The total mass values and uncertainties at interesting 
radii are derived from the 1000 simulated density and tem- 
perature profiles. The determination of the best-fit value 
and the la errors (at two sides) is mentioned in §4. One 
difference from the determination of entropy is that we in- 
clude only simulations that produce physically meaning- 
ful mass density profiles (mass density always larger than 
zero). The derived M500 — T500 relation with the tier 1+2 
groups is shown in Fig. 15. The BCES fits are listed in Ta- 
ble 6. We also included 14 clusters from V08 to constrain 
the M — T relation in a wider mass range. Our results 
show that the A/soo — Tsoo relation can be described by 
a single power law down to at least M500 = 2 x 10"'^^ 
h~^ M0. At T500 < 1 keV, more systems are needed 
to examine whether the relation steepens or not. Our 
-^500 — T500 relation is steeper than but still consistent 
with V08's (1.65±0.04vs. 1.53±0.08). Our slope is consis- 
tent with the Borgani et al. (2004) simulations with non- 
gravitational processes included (1.59±0.05), especially if 
two groups at T500 < 1 keV are excluded. The derived 
-^^500 ~ ^500 relation can also be compared with that by 
Arnaud et al. (2005), which has a slope of 1.71±0.09. Ar- 
naud et al. (2005) defined the system temperature as the 
overall spectroscopic temperature of the 0.1 r2oo — 0.5r2oo 
region, which should be close to Tiooo defined in this work 
(note rsoo ~ 0.66r2oo)- The data of the tier 1 and 2 groups 
give Tsoo/Tiooo = 0.96 on average. Then we find that the 
Chandra Afsoo — T relation constrained from 23 groups -|- 
14 V08 clusters is 18% - 3% higher than that by Arnaud et 
al. (2005) at 1 - 10 keV. We also notice that at 1 keV, the 
normalization of the Chandra Afsoo — T500 relation is 54% 
higher than that by Finoguenov et al. (2001). This is ex- 
pected from the generally higher density gradient around 
r5oo derived in this work (Fig. 16) than the typical val- 
ues in Finoguenov et al. (2001) (also see the Appendix of 
V06). 

As shown in Fig. 15, our M500 — Tsoo relation is offset 
from the true M500 — T500 relation in the NKV07 simu- 
lations, from 33% lower at 1 keV to 9% lower at 10 kcV. 
The agreement with the Borgani et al. (2004) simulations 
is better, from 18% lower at 1 keV to 6% lower at 10 
keV. Interestingly, the Msoo.hse — T500 relation from the 
NKV07 simulations is almost the same as the Chandra 
-^500 ^ 7500 relation from 23 groups + 14 V08 clusters. 
NKV07 attribute turbulence as the extra pressure to devi- 
ate the HSE mass from the true mass (also see Rasia et al. 
2004; Kay et al. 2004). Indeed for the ICM without a mag- 
netic field, the dynamic viscosity is roughly proportional to 
TjqIj. Thus, it may not be surprising that cool groups can 
develop stronger turbulence. However, the ICM is mag- 
netized and the real magnitude of the ICM turbulence is 
unknown. The NKV07 simulations only have numerical 
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viscosity that is small. Simulations with viscosity at dif- 
ferent strengths are required to better determine this bias 
term. High resolution X-ray spectra of the ICM may be 
ultimately required to constrain the turbulence pressure 
in the ICM. 

Kravtsov et al. (2006) suggested a new mass proxy, the 
Yx parameter (product of the gas temperature and the 
gas mass derived from the X-ray image, or Mgas, 500^500 
in this work), which in simulated clusters has a remark- 
ably low scatter of only 5%-7%, regardless of whether the 
clusters are relaxed or not. The agreement between sim- 
ulations and observations is also better for the M — Yx 
relation than the M - T relation (NKVG7). We exam- 
ined the M500 — Yx.5oo relation for 23 groups 4- 14 V08 
clusters (Fig. 15 and Table 7). Our results indicate 
that a single power law relation can fit the data very 
well. Our best-fit (0.57±0.01) is the same as the V08 
best-fit at Fx > 2 X IQI^Mq keV (0.57±0.G5), imply- 
ing the groups aligned well with clusters. Our best-fit 
is also consistent with the XMM-Newton result by Ar- 
naud et al. (2007) (10 clusters at Fx > lO^M© keV), 
on both the slope (0.548±0.027) and the normalization 
(within ~ 3%). Maughan et al. (2007) assembled 12 
clusters at z = 0.14 — 0.6 from their work and literature 
(Fx > 8 X IO^^Mq keV) and found that the slope of the 
-^500 — ^,500 relation is consistent with the fit to the V06 
clusters (o!564±0.009). Intrinsic scatter in both the M-T 
and Al — Yx relations are consistent with zero as the; mea- 
surement errors for groups are large. If we simply move the 
best-fit lines up or down to estimate the range of the scat- 
ter from the best-fit mass values, the scatter in the M — Yx 
relation is about the half that in the M — T relation. The 
slope is very close to the self-similar value of 0.6, especially 
if only tier 1 groups are included (0.588±0.012). Our best 
fits lie between the true mass and the HSE mass from the 
NKV07 simulations, but the offset is much smaller than 
that in the M500 — Tsoo relation. Thus, the Yx parame- 
ter appears to be a robust mass proxy down to at least 
2 xl0^3 1^-1 -y^i^i^ ^Yie derived M - T and M - Yx 

relations in this work, the slope of the Mgas,5oo ^Tsoo rela- 
tion is 1.89±0.05 (for tier 1+2 groups + ¥08 clusters), or 
1.86±0.06 (for tier 1 groups -|- V08 clusters). This value 
of the slope is consistent with the result by Mohr et al. 
(1999), 1.98±0.18 (90% confidence). 

7.2. Gas fraction 

We derived the enclosed gas fraction profile for each 
group. The enclosed gas fraction generally increases with 
radius and this trend continues to the outermost radius in 
our analysis, as generally found in V06. The enclosed gas 
fractions for groups at r25oo have a large scatter (Fig. 17). 
For groups with similar Tsoo, /gas, 2500 can be different by 
a factor as large as 2.5. Both the weighted average and the 
median of /gas, 2500 in our sample is 0.043, much smaller 
than the typical value of ^ 0.09 for V08 clusters. This 
mean can be compared with the average /gas, 2500 by G07 
(0.050±0.011), and we note that the 16 groups in G07 are 
on average brighter than our systems in a similar redshift 
range. The intrinsic scatter of /gas, 2500 is tightly correlated 
with the intrinsic scatter of i^o.iSrSOO) after the tempera- 
ture dependence of both variables are removed from their 
relations with temperature (Fig. 18). Groups with low 



/gas, 2500 have high i^o.isrsoo, relative to the average rela- 
tions. Thus, the large scatter of gas fraction within r25oo 
for groups is tightly correlated with the large scatter of en- 
tropy at 0.15 rsoo, and likely also the large scatter of X-ray 
luminosities for groups. Group properties (e.g., luminosity 
and central entropy) have large scatter because of the large 
scatter of the gas fraction around the center (e.g., r25oo)- 
Groups are on average fainter than what is expected from 
the self-similar L — T relation because groups are generally 
"gas-poor" within r25oo, compared with clusters. 

The enclosed gas fraction within r5oo is also derived for 
tier 1-1-2 groups (Fig. 19). We added 14 clusters from 
V06 and V08 to constrain the /gas — T500 relations. The 
/gas,500 ^Tsoo relation has a slope of ^ 0.16 - 0.22, depend- 
ing on whether both tier 1 and 2 are included (Fig. 19). 
We also give the /gas, 500 — Afsoo relation from the BCES 
orthogonal fit. For tier 1-1-2 groups -|- clusters: 



/gas,500 = (0.06 16 ± 0.0060) /if3^-^( 



1.5/ -^500 \0.135±0.030 



lOl3/lf3^M0 



For tier 1 groups -|- clusters: 



(8) 



/gas,500 = (0.0724 ±0.0078) ^^I'i ^^,^"1^/ '''^'-''' 

(9) 

We notice that the NKV07 simulations predict /gas, 500 oc 



rOoo52 (oj. mOoo"), from their best-fit M500 - T500 and 
Afsoo^Fx.soo relations. The gas fraction predicted in simu- 
lations depends on the modeling of cooling (e.g., Kravtsov 
et al. 2005) and is often tangled with the problem of pre- 
dicting the right stellar mass fraction in clusters. We also 
derived the enclosed gas fraction between r25oo and rsoo 
for 23 groups in tiers 1 and 2 (Fig. 19). Combined with the 
V06 and V08 results for clusters, the average /gas — 2500 
has little or no temperature dependence with an average 
value of ~ 0.12, although the measurement errors are not 
small, /gas, 2500-500 can also be derived as: 



/ 



gas, 2500-500 



1 



gas, 500 



-) 



(10) 



1 -a 

(a = M250o/Af500 = 5(?'2500/?'500)^) 

We use: /gas,25oo = 0.0347T|'o^°9 (from the BCES fit 
to 43 groups and 14 clusters. Fig. 17) and fgas 500 = 
0.0708r5"opo (from the BCES fit to 23 groups and 14 
clusters. Fig. 19). Combining our results on C500 
(§7.3) with V06's for clusters, roughly we have C500 = 
5.0(M5oo/10i3jy^o)"°-°®- Assuming an NFW profile, the 
^2500/^500 ratio can be well approximated as: ^2500/^500 ~ 
0.322 + 0.178 /5(1.523c5oo) at C500 = 1.3-5.0. Thus, com- 
bined with the M500 — Tkoo relation derived in this work, 
we can estimate the /gas, 2500-500 — Tsoo relation. Indeed 
as shown in Fig. 19, /gas, 2500-500 is nearly constant at 1 
- 10 kcV. The average /gas. 2500-500 is still ^ 27% lower 
than the universal baryon fraction (0.1669±0.0063 from 
Komatsu et al. 2008). However, one should be aware that 
the enclosed gas fraction still rises beyond rsoo, as gener- 
ally found in V06 and this work. We also notice that the 
observed /gas, 2500-500 is consistent with what was found 
in the simulations of Kravtsov et al. (2005) with cooling 
and star formation. We conclude that the low gas fraction 
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generally observed in groups is mainly due to the low gas 
fraction of groups within r25oo , or the generally weak abil- 
ity of the group gas to stay within r25oo • Beyond ^2500 , the 
groups are more regular and more similar to hot clusters, 
as also shown by the smaller scatter of their entropy values 
at r >r25oo- 

A natural question motivated by these results on en- 
tropy and gas fraction is: what is the fraction of the group 
luminosities within 0.15 rsoo or ^2500? Detailed work on 
the group Lx — T relation is beyond the scope of this pa- 
per and will be presented in a subsequent paper with an 
extended sample (including non-relaxed groups). Here we 
give the curves of the enclosed count fluxes for 17 groups 
that r5oo is reached by Chandra or (and) PSPC (Fig. 20). 
As the previous Lx~T relations only used the global spec- 
trum to convert the count rate to the group flux, Fig. 20 
can be regarded as the growth curves of the enclosed group 
luminosities. At 0.15 rsoo, the fraction ranges from 13% 
to 69%. At r25oo, the fraction ranges from 51% to 94%. 
Two extreme cases defining the boundaries are A1238 and 
AS 1101, also two systems with similar T500 but very dif- 
ferent /gas, 2500 (Table 3). A system with a bright cool 
core (like ASllOl) has a large fraction of the X-ray lumi- 
nosity within 0.15 rspo. Its system temperature without 
excluding the central core (e.g., 0.15 rsoo) therefore is bi- 
ased to a lower value. A system without a bright cool core 
(like A1238) has a small fraction of the X-ray luminosity 
within 0.15 r5oo, and its system temperature without ex- 
cluding the central core (e.g., 0.15 r5oo) is usually biased 
to a higher value (Fig. 3-5). All these factors contribute 
to the large scatter of the group Lx — T relation (e.g., Os- 
mond et al. 2004). Thus, the group Lx — T relation can 
be significantly contaminated by the large difference in the 
cores (e.g., within 0.15 rsoo). 

7.3. C500 

We also fitted the total mass density profile with the 
NFW profile and derived the concentration parameter, 
C500 = '^500 /fs , where Ts is the characteristic radius of the 
NFW profile. V06 used an inner radius of 0.05 rsoo, since 
the stellar mass of the cD is dominant in the center. The 
V06 sample is mainly composed of clusters. The groups 
in our sample have rsoo of 440 - 800 kpc, so the contri- 
bution of the stellar mass at 0.05 rsoo is still significant 
for low temperature systems (see e.g., G07). Thus, we use 
a fixed inner radius of 40 kpc. The outer radius is the 
outermost radius for the spectral analysis (rdot.spe in Ta- 
ble 2). We derived the total mass with 1 a uncertainties 
at radii corresponding to the boundaries of radial bins for 
spectral analysis between 40 kpc and rdct.spc (Fig. 3-5). 
The resulting mass density profile (at 4-10 radial points 
in this work) is fitted with an NFW profile and the un- 
certainty is estimated from 1000 Monte Carlo simulations. 
The results for 33 groups are present in Table 3 and are 
plotted with the system mass in Fig. 21. For the other 
10 groups, C500 is very poorly constrained. Our errors on 
C500 are larger than those in V06 and G07 as our errors on 
mass are larger. 

We compare our results with the best-fit c — M rela- 
tion from G07. G07 gave the best-fit Cvir — Afvir rela- 
tion. For their range of Cvir, C500 ~ 0.51cvir (for Cvir 

^ http:/ /www. physics. ucLedu/'bullock/CVIR/ 



= 10.35). For A = 101, we convert G07's best-fit re- 
lation to: C5oo(l + z) = 3.96(A/5oo/lOi'*M0)-"-226 (ad- 
justed to our cosmology). As shown in Fig. 21, at 
M500 > 4.5 X lO^^Afo, the G07 fit describes our results 
very well. But our results do not show significant mass 
dependence, so at 1/500 < 4.5 x 10^"^ Mq, our results are 
systematically below the G07 fit, although the errors are 
not small. However, the difference mainly comes from 
three groups (NGC 1550, NGC 533 and NGC 5129) for 
which G07 found — Al - 46 kpc (adjusted to our cos- 
mology), and our inner radius cut at 40 kpc prevents us 
from measuring such a small value of rg. In fact, the over- 
density radii of these three groups agree better between 
G07 and this work: rjsoo = 206 ± 2 kpc (G07) vs. 222±6 
from this work for NGC 1550, ri25o = 251 ± 2 kpc (G07) 
vs. 275±30 for NGC 533, ri25o = 217 ± 7 kpc (G07) vs. 
236±13 for NGC 5129. Excluding these three groups, our 
results agree well with G07's. Thus, the difference mainly 
hinges on the determination of rg that is related to the 
subtraction of stellar mass, while the results at large radii 
agree better. We also compare our results with the simu- 
lations of Bullock et al. (2001) and further workQ. As C500 
is sensitive to the halo formation time, smaller erg, SIm and 
tilt drive C500 smaller. We used the parameters derived in 
Komatsu et al. (2008) (see the caption of Fig. 21). As 
shown in Fig. 21, our results are generally consistent with 
the Bullock et al. simulations. Detailed discussions on the 
difference of the observed concentration parameter from 
the prediction can be found in Buote et al. (2007). We 
should point out that both our analysis and the V06 anal- 
ysis do not subtract the stellar mass and the X-ray gas 
mass, while G07 subtracted both components. G07 also 
included the XMM-Newton data for most of their groups. 
Future work on C500 may need very good measurement of 
the gas properties to > rsoo and careful modeling of the 
stellar and gas components (e.g., G07). The group dark 
matter mass profile may also need to be examined first to 
see whether a single NFW profile is the best fit. 

7.4. The baryon fraction and fossil groups in this sample 

We can estimate the enclosed baryon fraction from 
the cluster stellar fraction estimated before. Lin et al. 
(2003) gave: M,,5oo = 7.30 x lO"M0(rx/lkeV)i i69 
(adjusted to our cosmology) from the 2MASS data of 
nearby groups and clusters, which includes the stellar 
mass in cluster galaxies. With the M — T relation de- 
rived in this work from 23 groups and 14 V08 clus- 
ters, /*,500,Lin = 0.0263(Tx/lkeV)-0-4*i. Gonzalez et 
al. (2007) included the intracluster stellar mass and 
found: /*,5oo,Go„zaicz = 0.0380(Af50Q/10i4MQ)-0-64 (ad- 
justed to our cosmology). With our M50Q — T^oo relation, 
/*,500,Gonzaioz = 0.0864(T5oo/lkeV)-i o56.^ Adding the re- 
lation for /gas, 500 : the total baryon fraction from groups 
to clusters can be estimated. As shown in Fig. 19, there 
is substantial difference for groups. We also examined the 
stellar mass of the cD and its relation with other group 
properties. As shown in Fig. 22, the stellar mass of the 
group cD (which is proportional to its Kg band luminos- 
ity) is weakly correlated with the system mass. Low-mass 
cDs generally reside in low-mass groups. We also exam- 
ined the relation between Lks of the cD and /gas, 2500 but 
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found no correlation. However, in low-mass systems with 
low /gas, 2500) the Stellar mass of the cD can be comparable 
to the gas mass within r25oo- 

We also searched for fossil groups in this sample. Prom 
Jones ct al. (2003), fossil groups arc defined as a bound 
system of galaxies with the R band magnitude difference 
of the two brightest galaxies within half the virial radius 
larger than 2 mag. In this work, because we do not have 
homogeneous R band magnitudes for group galaxies in this 
sample, we used the 2MASS Kg band magnitude, which 
is a good indicator of the stellar mass. To ensure large 
and blue spirals are not left out, we also checked NED and 
HyperLeda to examine the B band magnitude difference. 
Jones et al. (2003) used r2oo as the virial radius, while 
we use the exact definition of the virial radius for our cos- 
mology, Tyir = (A ~ 101). Por the typical mass con- 
centration of groups in this sample, 0.5 rvir ~ r^oo- Thus, 
we examined the Kg and B band magnitude difference for 
group galaxies within r^ao- Six fossil groups arc selected: 
NGC 741, ESQ 306-017, RXJ 1159+5531, NGC 3402, 
ESQ 552-020 and ESQ 351-021. RXJ 1159+5531 and 
ESO 306-017 are known fossil groups (Vikhlinin et al. 
1999; Sun ct al. 2004). We note that NGC 1132 was 
considered a fossil group (Mulchaey & Zabludoff, 1999). 
However, NGC 1126 is 8.4' from NGC 1132 with a veloc- 
ity difference of 438 km/s (r5oo = 16.2' or 440 kpc for the 
NGC 1132 group). The 2MASS magnitude difference 
is 1.49 mag and the B magnitude difference is 2.07 mag. 
Thus, we do not consider NGC 1132 a fossil group based 
on the Jones et al. (2003) definition. 

Although the sample is small and not representative, we 
examined whether these fossil groups populate a different 
position of the scaling relations than non-fossil groups. No 
significant difference in K — T, /gas — T and M — T rela- 
tions is found. As shown in Fig. 22, these fossil groups 
indeed have massive cDs. The C500 of these fossil groups 
(four listed in Table 3, others poorly constrained) are ^ 
2.4 - 4.3, on average smaller than the average of this sam- 
ple, which differs from the claim by Khosroshahi et al. 
(2007) that fossil groups have higher mass concentration 
than non-fossil systems. 

8. AGN HEATING AND RADIO GALAXIES 

8.1. Sign of heating in entropy profiles 

Besides prc-heating (from SNe and AGN) and cooling, 
impulsive heating from the central AGN is often required 
to explain the observed scaling relations like L — T and 
K — T (e.g., Lapi et al. 2005). Strong shocks driven by the 
central AGN may boost the ICM entropy and create an en- 
tropy bump. This transient anomaly in the entropy profile 
may be detected in a large sample. We have searched for 
such entropy features in our sample and find two promis- 
ing cases: UGC 2755 and 3C 449 (Pig. 23). There are also 
three groups with a significant break observed in their sur- 
face brightness profiles (3C442A, IC1262 and A2462) that 
may be related to AGN heating. We briefly discuss them in 
the Appendix. Both UGC 2755 and 3C 449 host a strong 
PR I radio source with two-sided radio lobes. UGC 2755 
has a central corona with a radius of ^ 3 kpc (see Sun 
et al. 2007 for the connection of thermal coronae with 
strong radio sources). Prom 10 kpc to 80 kpc in radius, 
the surface brightness profile is very flat. Then there is 



a sharp break at ^ 90 kpc. UGC 2755's radio lobes ex- 
tend to ^ 100 kpc from the nucleus in the NVSS image, 
which may naturally explain the entropy bump within ~ 
90 kpc. 3C449 also has a central corona with a radius 
of ^ 3 kpc. Its entropy bmnp is at ^ 40 kpc - 100 kpc 
and less significant as that in UGC 2755. 3C449's radio 
outflow is spectacular and can be traced to at least 1.6 
Mpc in radius from the NVSS image (compared to its r5oo 
of 433 kpc estimated from r — T relation). Nevertheless, 
the brighter inner part of the radio jets/lobes ends at ^ 
100 kpc from the nucleus. Thus, the entropy bump we 
observed in 3C 449 may represent the most recent heating 
event. 

Strong shocks are required to effectively boost entropy. 
Using the standard Rankine-Hugoniot conditions, the en- 
tropy increase after shocks is a sensitive function of the 
shock Mach number, and weak shocks have little eSect on 
amplifying entropy. A shock with a Mach number of 1.2 
only increases entropy by 0.44%. A Mach 2 shock increases 
entropy by 20%, while Mach 3.3 and 4.4 shocks increase en- 
tropy by 100% and 200% respectively. The entropy bump 
in 3C449 is only at the level of 20%, which can be pro- 
duced by a single Mach 2 shock. The entropy boost in 
UGC 2755 is 1.5 - 2 times, which can be produced by a sin- 
gle Mach 3 shock. The adiabatic sound speed in groups is 
not high, 540 - 630 km/s in 1.1 - 1.5 keV ICM for these two 
groups. In these low density groups, the ambient pressure 
is much lower than that in the dense cores of hot clusters. 
Shock deceleration may also be slower. The velocity of the 
outflow-driven shock is ~ fp{Pki-a/ pr'^Y^^ , where p is the 
ICM density and fp is a structure factor of order unity 
that depends on the outflow geometry and the preshock 
density profile (e.g., Ostriker & McKee 1988). rig ~ 10~^ 
cm~^ around the entropy bumps and Pkin ~ 10'^'' ergs 
s~^ from the radio luminosities of two radio AGN and the 
relation derived by Birzan et al. (2004). The estimated 
velocity is then ~ 1300/p km/s, which is comparable to 
the requirement of the entropy boost. Therefore, radio 
outflows in these two groups are capable of driving Mach 
2-3 shocks to produce the observed entropy bumps. 

8.2. Strong central radio sources in this sample 

There are ten groups hosting a central radio source that 
is more luminous than Z/i.4ghz = 10^^ W Hz~^: 3C 31, 
3C 449, UGC 2755, 3C 442A, A160, A2717, ASllOl, 
A3880, A1238 and A2462. What is the typical X-ray gas 
environment around these radio sources? Six of them lack 
large cool cores (e.g., <; 30 kpc radius), 3C31, 3C 449, 
UGC 2755, A160, A1238 and A2462. However, all of them 
host a central corona with a radius of ^ 3 - 8 kpc, typi- 
cal for massive cluster and group galaxies as discussed in 
Sun et al. (2007). This component is reflected in their 
temperature profile, except for A1238 as it is faint. The 
spectrum of A1238's central source can be described by a 
^ 0.8 keV thermal component. Its X-ray luminosity and 
Kg band luminosity fall on the typical region for coronae 
and its properties are similar to ESO 137-006 in A3627 
(a nearby bright corona. Sun ct al. 2007). Four other 
groups (3C442A, A2717, A3880 and ASllOl) host larger 
cool cores with a radius of ^ 30 kpc or larger. The cool 
cores in 3C442A and A2717 are clearly disrupted, likely 
by the radio sources. Thus, all these strong radio sources 
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have low-entropy ICM (< 30 keV cm^) at the center. 

9. SYSTEMATIC ERRORS 

We follow Humphrey et al. (2006), GOT and V08 to 
discuss the systematic error budgets in our results. The 
uncertainties of the local X-ray background arc the main 
error budget at large radii. As shown in §3.2, we have 
included a conservative estimate of the background uncer- 
tainties into the errors of temperatures and densities. This 
is the primary reason that our results at large radii have 
larger errors than those of V06 for the same systems. 

We used the LAB survey data (Kalberla et al. 2005) 
for the galactic hydrogen column density (Table 2) and 
examined the Chandra spectra for excess absorption. The 
LAB column density is on average ~ 6% lower than the 
column density from Dickey & Lockman (1990). Out of 
43 groups we studied, 19 show significant excess absorp- 
tion. We can compare this fraction to that of V05, 6 out 
of 13 with excess absorption. For the five groups that 
were studied in V05 and this work, both works find ex- 
cess absorption for the same three groups. On the other 
hand, GOT used the galactic hydrogen column density from 
Dickey & Lockman (1990) for all groups. The incidence 
of excess absorption in our work increases with the galac- 
tic hydrogen column density. At A^2icm > 4 x 10^° cm~^, 
eight of eleven groups (most of them at z < 0.03) show 
excess absorption. This trend is qualitatively consistent 
with the result by Arabadjis & Bregman (1999). We also 
examined the effects of a conservative Nn uncertainty of 
±2 X 10^" cm~2 (see e.g., V08) on our results. For T < 
1.6 keV gas, the determination of temperature is little af- 
fected by absorption. An TVh change of ±2 x 10^" affects 
the temperature by ~ =F 1-4% and the density by ~ ± 
6.1%. The subsequent changes on K,M,Yx and /gas at 
an overdensity radius are ~ =f6.3%, =f2.1%, ±3.9% and 
±T.6% respectively, assuming an entropy slope of O.T (note 
that the overdensity radius also depends on mass). For 
T > 1.6 keV gas, an A^h change of ±2 x 10^° affects the 
temperature by ~ =F 5.3% and the density by ± 3.T%. 
The subsequent changes on K, M, Yx and /gas at an over- 
density radius are ~ t5-9%, t8.0%, t4.4% and ±9.T% 
respectively. At small radii, uncertainties at these levels 
are not important as they are smaller than the intrinsic 
scatter (e.g., ii'o.isrsoo and /gas, 2500)- At large radn, the 
statistical errors (including uncertainties from the local X- 
ray background) overwhelm. Nevertheless, the systematic 
errors from A^h should be kept in mind. 

We use the deprojection algorithm derived by Vikhlinin 
(2006). The form of the 3D temperature profile is the same 
as that used in V06. The robustness of this deprojection 
algorithm has been presented in Vikhlinin (2006), with un- 
certainties of ^ 0.05 keV. Nagai et al. (200Ta) presented 
mock Chandra analysis of cluster simulations, using the 
deprojection algorithm proposed by Vikhlinin (2006) and 
the form of the temperature profile used in V06. The best- 
fit 3D temperature profile is well consistent with the true 
temperature profile in simulations, with residuals in a sim- 
ilar level as shown in Vikhlinin (2006). This test further 
validates the approach used in the present work. We have 
also performed a test on a sub-sample of 6 groups with 
best-quality data. The traditional onion-peeling method 
(e.g., McLaughlin 1999; GOT) was used. The resulting de- 



projected temperature profiles often show small magnitude 
of oscillating. If we use the form of equ (2) to fit them, 
the best-fits are consistent with those shown in Fig. 3-5 
within ~ 0.05 keV on average. The uncertainties on tem- 
peratures range from 2% - 6%, while the uncertainties on 
densities are 0.5% - 3% (both directions). Thus, this sys- 
tematic error (larger for < 1 keV groups) are smaller than 
the systematic error related to the uncertainty of A^h- 

There are other systematic errors, like the choice of the 
plasma spectral codes, and departures from spherical sym- 
metry in the group gas. Those are either small in magni- 
tude or have little impact on the scaling relations, as dis- 
cussed in Humphrey et al. (2006), GOT and V08. As mock 
data from simulations arc starting to be analyzed in the 
same way as the observational data (Nagai et al. 200Ta; 
Rasia et al. 2008), these factors are becoming better con- 
trolled. 

10. SUMMARY AND CONCLUSIONS 

We present an analysis of 43 galaxy groups with Chan- 
dra observations. With inclusion of many faint systems 
(e.g., the ones hosting strong radio sources), our sample is 
not much biased to X-ray luminous groups, as shown by 
the wide range of ICM entropy values around the group 
center (§6). We used the ACIS stowed background and 
modeled the local CXB for each group. Uncertainties of 
local background arc folded into the derived temperature 
and density profiles. The projected temperature profile 
and the surface brightness profile are modeled with sophis- 
ticated models, which have enough freedom to describe the 
data from the core to the outskirts. The 3D abundance 
profile is also derived. The 3D temperature and density 
profiles are constrained through iterative fitting. The un- 
certainties are estimated from 1000 Monte Carlo simula- 
tions. Gas properties are derived to at least r25oo for all 
43 groups. For 11 groups, we can derive gas properties 
robustly to rsoo. For another 12 groups, gas properties 
can be robustly derived to ^ riooo, so we extrapolate the 
results to rsoo- The main results of this paper are: 

1) We present the M500 — T500 and M500 — Yx.^oo rela- 
tions in M500 = IO^^H'^Mq - lO^^/i"^M0, combined with 
the V08 results on 14 T500 > 3.T keV clusters. Both re- 
lations are well behaved at the low-mass end and can be 
well fitted with a single power law (§T.l and Fig. 15). The 
-^500 ~ ^.500 relation indeed has a smaller scatter than 
the M500 — Tsoo relation (about half). The M500 — T500 
relation from observations is still offset from simulations 
(e.g., NKVOT). Interestingly, the Chandra M500 — Tsoo re- 
lation is very close to the Mboo.hse — T'500 relation in the 
NKVOT simulations. Although it is tempting to attribute 
the difference to a mass bias, better understanding of the 
ICM viscosity is required. 

2) The group gas fraction within r25oo is on average 
much smaller than that of clusters (e.g., ~ 0.043 for ~ 1.5 
kcV groups vs. ^ 0.09 for ^ 6 keV clusters, §T.2 and Fig. 
IT), which is consistent with GOT's results. The group gas 
fraction within r25oo also has a large scatter, spanning a 
factor of ^ 2 at any fixed temperature. On the other hand, 
the gas fraction measured between r25oo and rsoo has no 
temperature dependence with an average value of ~ 0.12 
(Fig. 19). Thus, the generally low gas fraction in groups 
is due to the general low gas fraction within r25oo- 
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3) We derived the K — T relations at 30 kpc, 0.15 rsoo, 
?'2500, J'lsoo, r-iooo and rsoo (§6.1, Fig. 11). The large scat- 
ter of the entropy values at 30 kpc and 0.15 rsoo reflect 
the wide luminosity range of groups in this sample. The 
K — T relation is significantly tighter beyond r25oo and 
the intrinsic scatter of entropy is the same at 10% - 11% 
from r25oo to rsoo- Thus, the group properties are more 
regular from r25oo outward, in line with the gas fraction 
results. With 14 clusters from V06 and V08 included, wo 
also present K — T relations in the full temperature range 
(Fig. 12 and Table 5). At rsoo, the slope of the K-T rela- 
tion is consistent with the value from self-similar relation 
(1.0). 

4) The ratios of the observed entropy values to the base- 
line values (from adiabatic simulations) decrease with ra- 
dius (§6.2 and Fig. 13). At r2soo, the ratio ranges from 
~ 1.8 - 3.5, with a weighted mean of 2.2. The weighted 
mean decreases to ^ 1.8 at riooo and ^ 1.6 at rsoo. The 
still significant entropy excess at rsoo in groups may re- 
quire a diffuse way to distribute heat (e.g., AGN heating, 
see §8.1), but it may also be understood with smoother ac- 
cretion and the mass bias that may be especially large in 
groups (§6.2). The entropy excess at rsoo is also detected 
for 14 clusters from V06 and V08 (~ 35%, §6.2). The dif- 
ference in the entropy excess at rsoo between groups and 
clusters (~ 17%) is not as large as previously claimed from 
the ROSAT and ASCA data. 

5) The entropy slopes arc determined at 30 kpc - 0.15 

rsoo, 0.15 rsoo - ''asoo, ?'250o - ?'i500 and risoo - ''det.spe (§6.3 
and Fig. 14) . The slopes are all shallower than 1.1 be- 
yond 0.15 rsoo- Scatter is large but the average slope is ~ 
0.7 beyond 0.15 rsoo. 

6) The group temperature profiles are similar at > 
O.lSrsoo, despite large scatter within 0.15 rsoo (§5 and 
Fig. 8). The average slope is consistent with the "Uni- 
versal temperature profile" (Markevitch et al. 1998; De 
Grandi & Molendi 2002; Loken et al. 2002; V05; Pratt et 
al. 2007; G07; Rasmussen & Ponman 2007; LM08) but 



there is still scatter. The group temperature profiles also 
appear more peaky than those of clusters (Fig. 9). 

7) We also derived the concentration parameter (csoo) 
for 33 groups (§7.3). Our results are generally consistent 
with the Bullock et al. (2001) simulations under the cur- 
rent WMAP5-FSN-hBA0 cosmology. 

8) We selected six fossil groups in this sample and four 
are new. The X-ray gas properties of these fossil groups 
have no significant difference from non-fossil groups (§7.4) 
in scaling relations. 

9) We found two groups with substantial entropy bumps 
(§8.1), which may indicate a recent strong heating episode. 
Both host strong radio galaxies at the center and we esti- 
mate that the radio AGN is capable of driving shocks to 
boost entropy to the observed level. 

The emerging picture of groups from this work is that 
the main difference between groups and hotter clusters is 
the general weak ability of the group gas to stay within 
''2500 , which explains most of deviation of the group prop- 
erties from the self-similar relations (e.g., entropy and lu- 
minosity). The group properties within r25oo have large 
scatter, but may all be related to variations in the level 
of the enclosed gas fraction within r25oo- Beyond r25oo, 
groups are more regular and more like clusters, making 
them promising tools for cosmology, as shown by the well- 
behaved M — T and M — Yx relations derived in this work. 

This work would be impossible without the rich data 
sets in the Chandra archive. We would like to thank 
all Pis of the Chandra observations for their original ef- 
forts. We thank M. Markevitch for his help on the ACIS 
stowed background. We thank the referee, D. Buote, for 
helpful and prompt comments. We thank helpful discus- 
sions with M. Balogh, S. Borgani, D. Nagai and G. Pratt. 
The financial support for this work was provided by the 
NASA LTSA grant NNG-05GD82G. We made use of the 
NASA/IPAC Extragalactic Database and the HyperLeda 
database. 



APPENDIX 
Components of the Chandra background 

The Chandra background has been extensively discussed before (e.g., Markevitch et al. 2003; Wargelin et al. 2004; 
V05; HM06; Humphrey & Buote 2006). We present here a brief summary, focusing in particular on the cosmic X-ray 
background (CXB). There are two basic components in the quiescent Chandra background, particles and photons (or 
CXB). The dominant background component is the charged particle background (PB), which is not vignetted. The 
spectrum of the Chandra PB has been remarkably stable since 2000 (Vikhlinin et al. 2005, V05 hereafter; Hickox & 
Markevitch 2006, HM06 hereafter), although the absolute flux changes with time and is related to the solar cycle. The 
only exception so far is for BI data after the middle of 2005, which is discussed in the next section. 

The cosmic hard X-ray background is considered to be composed of unresolved X-ray point sources, mostly AGN. This 
component can be described by an absorbed power-law with a photon index of ~ 1.5 (HM06). Its flux depends on the 
level of point source excision or the limiting flux for point sources. HM06 analyzed the two deepest Chandra fields at that 
time, Chandra deep field north (a combined clean exposure of 1.01 Ms) and south (a combined clean exposure of 0.57 
Ms). The unresolved hard X-ray background flux density is (3.4±1.7)xl0~^^ ergs s~^ cm~^ deg~^ in the 2-8 keV band, 
which represents the lower limit of the hard CXB flux in shorter Chandra observations. The X-ray logN — logS relation 
has been well studied allowing us to predict the average unresolved cosmic hard X-ray background below the point source 
limiting flux. Kim et al. (2007; K07 hereafter) presented the average relation between point source limiting flux and 
expected hard CXB flux density (Fig. 19 of K07). For the limiting flux for point sources in the outermost bins of groups 
in our sample, we expect a resolved fraction of 25% - 65% in the 2-8 keV band, which corresponds to a flux density of 
6.1 - 13.1 xlO~^^ ergs cm~^ s~^ deg~^ for the unresolved hard CXB background. However, one should be aware that the 
K07 relation is just an average. About 30% variation can be expected over the angular scale we study (HM06). 
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The soft X-ray background is composed of several components, Galactic, local bubble, geocoronal and heliospheric 
emission (e.g., Wargelin et al. 2004; Snowden et al. 2004; Koutroumpa et al. 2007). The latter two components are 
primarily from the solar wind charge exchange (SWCX) emission and are time variable, with a contribution to the O 
VII and O VIII lines as much as the Galactic component (Koutroumpa et al. 2007). The strength of the soft X-ray 
background varies with the sky position (as shown in the R45, or 3/4 keV ROSAT all sky survey map). There has been 
lots of work done to quantify its spectral properties. With the Chandra data, Markevitch et al. (2003) and HM06 have 
shown that, the soft X-ray background beyond the regions with strong RASS R45 flux (e.g., the North polar spur), can 
be well described by a single unabsorbed thermal component with a temperature of ^ 0.2 keV. Its spectrum typically 
shows a broad line hump around 0.6 keV, mainly from the 0.57 keV O VII and 0.65 keV O VIII lines. The soft X-ray 
background has also been studied with the Suzaku data, which have the higher spectral resolution. Fujimoto et al. (2007) 
analyzed the Suzaku data of the North Ecliptic Pole region (R45 = 140 xlO~^ counts s~^ arcmin"^) and found that 
the soft X-ray component at the non-flare period has a temperature of 0.18 keV with over-solar abundances. Miller et 
al. (2008) analyzed the Suzaku data of the brightest region of the North Polar Spur (NPS, R45 = 748 xlO"^ counts 
s~^ arcmin"^) and found that the NPS thermal component has a temperature of ^ 0.28 keV with generally sub-solar 
abundances, besides the assumed 0.1 keV local bubble and galactic halo components with solar abundance. Snowden et 
al. (2008) adopted a soft CXB model composed of an unabsorbed ~ 0.1 keV component for the local hot bubble, an 
absorbed ^ 0.1 keV component for the cooler Galactic halo emission and an absorbed ~ 0.25 keV component for the 
hotter Galactic halo emission. All abundances are fixed at solar. Thus, we have enough knowledge to model the soft 
X-ray background. 

"Blank sky background" and stowed background 

Much previous Chandra work on clusters used the blank-sky background data0, which are good enough in high surface 
brightness regions. However, the averaged soft X-ray component in the standard blank-sky background data is very likely 
different from the actual soft X-ray foreground in any particular field, as pointed out previously by e.g., Markevitch et 
al. (2003), V05, Humphrey & Buote (2006) and G07. Taking the example of the blank sky background data in period 
D, the ACIS-I file combines 29 pointings from Jan., 2001 to Nov. 2004, at regions where R45=90 - 150. The exposures 
range from 27 ks to 165 ks with a median of ^ 66 ks. The ACIS-S file combines 12 pointings from Nov. 2001 to Oct. 
2003, at regions where R45=90 - 150. The exposures range from 20 ks to 114 ks with a median of ~ 30 ks. In 5 of 12 
pointings, there are no SI data. The ACIS response has been changed significantly from 2001 to 2004 (especially because 
of the contamination on the optical blocking filter). Thus, it is often inadequate to only use the blank-sky background 
to constrain the ICM properties at low surface brightness regions. Moreover, the PB and the CXB are not separated 
in the blank-sky background data. When the blank-sky background is scaled to account for the flux change of the PB, 
the CXB in the blank-sky background data is unphysically scaled. While the PB rate did not change much between 
the spring of 2000 and the spring of 2004, it has been significantly increasing ever since. In 2006, the PB rate was on 
average 50% higher that the average value between the spring of 2000 and the spring of 2004 (Fig. 6.24 from Chandra's 
Proposers' Observatory Guide v.Sfl). Thus, any analysis for data taken after the middle of 2004 involves a large scaling 
of the PB, often resulting significant over-subtraction of the CXB. Therefore, a second correction besides subtracting the 
scaled blank-sky background is required (e.g., Vikhlinin et al. 2005). This "double subtraction" is often efficient but 
requires the presence of source-free regions in the Chandra field, which is not true for many nearby groups in our sample. 
In this work, we utilize the newly available ACIS stowed background data to subtract the PB component. The CXB is 
modeled and the uncertainties are folded into the final error budgets. 

Since Sep., 2002, ACIS observations have been carried out twice a year in the stowed position, shielded from the sky 
by the science instrument module structure and away from the on-board calibration source. By the end of May, 2007, 
415 ksec data had been collected. Background flares have never been observed in the stowed data. The comparison 
with the dark moon observations indicates that the stowed background is the same as the quiescent PB collected by the 
CCDs in the normal focal position (Markevitch et al. 2003; the CXC calibration website ^ .). HM06 also used the stowed 
background (236 ks at the time of their work) to carry out absolute measurement of the unresolved CXB. They show 
that between Jan. 30, 2000 and Sep., 2002 when the stowed background data are not available, the spectral slope of 
the PB is the same. Thus, we can apply the stowed background to early data. The stowed background allows us to 
separate the non- vignetted PB from the vignetted CXB. The main reason for our preference for the stowed background 
over the blank-sky background is that we have better control of the local background for nearby groups where sources fill 
the whole Chandra field. The "double subtraction" method with the blank-sky background cannot be applied for these 
nearby groups, as there is no region that is free of group emission. The derived local X-ray background based on stowed 
background also has a clear physical meaning and can be compared between Chandra observations with very different 
PB fiuxes or observations with other telescopes like l^MM-Ncwton. Another subtle advantage of the stowed background 
is related to the telemetry limit of Chandra, especially for the VFAINT mode data. The blank-sky background data of 
each specific ACIS CCD may come from different combination of observations, especially for ACIS-S (e.g., SI vs. S2-I-I3). 
Therefore, the residual or decremental background on the SI CCD is in principle different from that on FI chips, which 
complicate the analysis. Similarly, even if S3 is turned on for ACIS-I observations, it cannot be used for local background 
study as the blank-sky background of the S3 chip (when the aimpoint is on ACIS-I) is only a subset of the ACIS-I 

* http:/ /cxc. harvard.edu/contrib/maxim/acisbg/ 
^ http:/ /asc. harvard.edu/proposer/POG/ 
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blank-sky background data. On the contrary, with the stowed background, we are analyzing the absolute CXB in the 
interested field so data on different CCDs can be fitted jointly to make better constraints. Therefore, we used the stowed 
background to subtract the PB in our observations. 

The spectra of the Chandra PB have been very stable (e.g., HM06). However, a small change of the spectral shape on 
the BI CCDs (S3 and SI) has been identified from aroimd the middle of 2005, while the spectra of the PB on the FI chips 
still keep the same ^. The change appears abrupt around the middle of 2005, while the spectral shape of BI data remain 
the same from that time to at least the middle of 2007 (private communication with Maxim Markevitch). The stowed 
background data have been broken into 2 periods, one with 235 ksec total exposures from five observations of Sep. 3, 2002 
to June, 10, 2005 (period D), the other with 180 ksec total exposure from four observations of Nov., 13, 2005 to May, 28, 
2007 (period E). We emphasize that the notation adopted here is only for the purpose of this work. We examined the 
spectral difi'erence between these two periods, after matching their fluxes in the 9.5 - 12 keV band. The flux of the period 
D background is always a little lower than that of the period E background, after re-scaling. The biggest difference is 
seen on the S3 chip, with a 6.0% difference in the 0.35-7.0 keV band (Fig. 24). The residual emission is very flat (note 
it is non- vignetted) and can be removed by increasing the D period background by 5.7%. The difference is smaller on 
the SI chip, 2.4% in the 0.35-5.5 keV band (note that the SI PB increases rapidly at > 5.5 keV, Fig. 6.21 of Chandra^s 
Proposers' Observatory Guide v. 9). The difference is consistent with zero in FI chips, ~ 1% in the 0.5-7.0 keV band, 
excluding the Au line in the 2.0-2.3 keV band. Thus, we can apply the total stowed background (415 ks exposures) to the 
FI data. For BI data, we use the stowed background in their corresponding periods. There are only two groups in our 
sample with BI data taken between June, 10, 2005 and Nov., 13, 2005 (NGC 1550 and NGC 5098). Both were observed 
after Oct. 22, 2005 and both have earlier FI data. We used the period E stowed background for the BI data of both 
groups. In this work, we also take a larger uncertainty on the normalization of period E PB for BI chips (5% compared 
to 3% for FI data and the period D BI data). 

NOTES ON SOME GROUPS 

In this section, we present notes on some groups, mainly on the comparison with previous work on the gas properties 
at large radii. Thus, the cited references are usually not complete for each group as the detailed dynamical and thermal 
structures of the group cores are beyond the scope of this work. 

NGC 1550 was examined by Sun et al. (2003), with two ACIS-I observations. Now with two additional longer 
ACIS-S exposures in the offset positions, the gas properties in this system can be constrained much better. NGC 1550's 
temperature profile is among the best determined for 1 keV groups, with the good Chandra coverage. We can compare 
our results with those from G07 who analyzed two short ACIS-I observations and an XMM-Newton observation. The 
temperature profiles agree well although we constrain the temperature to larger radii from the ACIS-S observations. Our 
rsoo and M500 are 10% and 50% higher than those derived by G07. Our C500 (4.931^;^^) is smaller than G07's (9.0±0.6). 

NGC 3402 has been studied by V05 and V06. Although our temperature profile is consistent with V05, the slope of 
the decline is smaller. The IRAS 100/im map shows the presence of Galactic cirrus around the group. We indeed derived 
a higher absorption column, l.lzbO.l x 10^^ cm~^, than the Galactic value (4.0x10^° cm~^ from LAB). This value is 
smaller than that derived in V05, 1.55±0.1 x 10^^ cm~^. However, the absorption difference has little effect on the derived 
gas temperature. As the gas temperature is mainly determined by blended line centroid, gas temperature remains almost 
the same, with the higher A^h in V05. Because of higher temperatures derived at large radii, our derived M2500 is larger 
than that from V06 and /gas, 2500 is smaller than in V06. 

Abell 262 is a nearby luminous system in which the X-ray emission can be traced to over 800 kpc in the 7.6 ks 
PSPC data. A262 was included in V06 and G07 samples. We included a new deep Chandra exposure (110 ks) taken in 
2006 in our analysis, while G07 also analyzed an XMM- Aewfon observation. Our results of 7'25oo , -^^2500 and /ga,5,2500 are 
consistent with those in V06 and G07. Our C500 (3.48I045) is well consistent with V06's result and is close to G07's result 
(4.5±0.4), but smaller than the result by Piffaretti et ai. (2005) (5.8±1.2). 

NGC 383 hosts a bright FRI radio source 3C31 . There is a backgromid cluster centered on 2MASX J01065891+3209285 
at z=0.1116. We derived the surface brightness profile centered on the backgroimd cluster and also analyzed a short 
exposure (ObsID 3555, 5.1 ks) centered on the background cluster. The cluster emission is detected to ~ 4' radius. In 
the analysis for NGC 383, we excluded the region within 6.5' of 2MASX J01065891-I-3209285, which is a bit larger than 
its r5oo (5.9' for kT = 2.3 keV). 

3C 449 is located at a Galactic latitude of -16 deg and at the outskirts of a bright IRAS 100 nm feature across several 
degrees, which should explain the enhanced absorption. It is one of the two groups with an entropy bump detected (§8.1). 

NGC 533 was studied by Piffaretti et al. (2005) (XMM-Newton) and G07 {Chandra + XMM-Newton). The explored 
radial range in spectral analysis is similar in all three work (up to 240 - 260 kpc). Our temperature profile agrees well 
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with that derived by G07. Our C500 (4.58l2 34) is stiH consistent with the results by Piffaretti et al (2005) (8.6±0.7) and 
G07 (9.0±0.7). 

MKW4 is bright enough that the group emission can be traced to the very edge of the Chandra field. We are able to 
separate the group emission and the soft background emission on the SI spectrum because of the prominent ~ 0.6 keV 
hump in the soft background emission and the iron L hump in the group emission. The same Chandra data had been 
analyzed by V05 and V06. Our temperature profile agrees well with V05's, as well as properties at r25oo (mass and gas 
fraction) with V06's. However, our results at rsoo differ from VG6's. The difference should lie on the modeling of the 
density profile. V06 derived a very steep density profile at large radii of MKW4 (/3eff,500 = 0.92, Table 2 of V06), while 
we derived a value of ^ 0.6, more similar to A262 and A1991 (Table 2 of V06). We notice that Vikhlinin et al. (1999) 
derived Pouter = 0.67±0.06 for MKW4, using the same PSPC data. G07 used the same Chandra data and also analyzed 
an XMM- Newton observation of MKW4. Our rsoo and M500 are consistent with G07's values, but our C500 lies between 
the values of V06 and G07. G07 assumed an NFW profile for the dark matter halo and only derived gas properties to 
322 kpc. The properties at rsoo thus rely on the assumption of NFW profile and extrapolation. On the other hand, 
the Chandra data at the outermost bin only covers 10% of the area in that annulus. Better constraints on the properties 
of this nearby system require more coverage at large radii. 

NGC 5129 is at the edge of the NPS so the local R45 value is very high. We indeed found a high local soft X-ray excess 
(Table 2). It is also near an extended feature on the IRAS 100 /im map, which should explain the enhanced absorption. 
G07 presented results based on an XMM-Newton observation. Our temperature profile is consistent with that from G07. 
However, our C500 is smaller than G07's (3.43ti;22 vs. 7.7±1.3). 

UGC 2755 is one of the faintest and most gas poor systems in our sample (/gas, 2500 = 0.030 ± 0.005). It is one of the 
two groups with an entropy bump detected (§8.1). 

NGC 4325 has a luminous cooling core. NGC 4325 was studied by G07, who also analyzed an XMM-Newton obser- 
vation. The explored radial range in spectral analysis is similar in both work (228 kpc vs. 232 kpc in our work). Our 
temperature profile agrees well with that derived by G07. Our C500, M2500 and /gas, 2500 are consistent with G07's within 
1 a errors. 

3C 442A has the most peaky temperature profile in this sample. The sharp reconstructed temperature peak at 75 - 
150 kpc, and the steepening of the surface brightness at ~ 120 kpc, may best be explained by the second most luminous 

radio source in our sample (only after Abell 2462). From the NVSS image, the two radio lobes of 3C 442A extend to ^ 5' 
(or 152 kpc) in radius. Shock heating by the radio source may explain the high temperature peak. Because of this high 
temperature peak, the total mass density profile is not physically meaningful within the central ^160 kpc, which casts 
doubt on the assumption of hydrostatic equilibrium within the central 160 kpc. Therefore, C500 cannot be constrained 
reliably. 

ESO 552-020 was also studied by G07. Our temperature profile is consistent with G07's. Our C500 is also consistent 
with G07's. 

IC 1262 has rich substructure within its core. However, beyond the central 60 kpc radius, it appears symmetrical and 
relaxed. Its surface brightness profile shows a sharp break at ~ 200 kpc. Unlike 3C 442A and A2462, its radio source is 
faint. However, this may be caused by past AGN activity, which was also suggested to explain the rich structures in the 
group core by Trinchieri et al. (2007). The sharp temperature decline beyond 200 kpc also supports this scenario, as the 
regions between 80 and 200 kpc radius may have been recently heated (Fig. 4). The adiabatic sound speed in ~ 1.9 keV 
gas is 700 km/s. Thus, a Mach 1.5 shock will travel to the current position in ~ 200 Myr, which is consistent with the 
typical duty cycle of radio AGN. 

ESO 306-017 has been studied before (Sun et al. 2004; G07). There is an adjacent Chandra pointing targeted at the 
z = 0.64 cluster RDCS J0542-4100 (ObsID 914) that we used to constrain the local soft CXB. Our temperature profile 
covers a wider radial range than G07's. In the overlapping region, our temperature profile is consistent with G07's. Our 
results on r^oo and M500 are also consistent with G07's. 

NGC 5098 was also studied by G07. The IRAS 100/im map shows the presence of Galactic cirrus around the group 
and we indeed find extra absorption (Table 2). There is a second group in the field (Mahdavi et al. 2005) and the region 
around it is excluded in our analysis. Our temperature profile and results on rsoo, M500 and C500 are consistent with 
G07's. 
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UGC 842 was also studied by G07. Our temperature profile and results on rsoo, M500 and C500 are consistent with 
G07's. 

A2717 was studied by Pratt & Arnaud (2005), GOT and Snowden et al. (2008) with the XMM-Newton data. Our 
temperature profile is consistent with the profiles derived by G07 and Snowden et al. (2008). Our derived rsoo and M500 
are consistent with those derived by G07, while our C500 (2.15^q |2) close to G07's value, 3.0±0.2 and the result by 
Pratt & Arnaud (2005), 2.8±0.2 (converted from their C200 assuming an NEW profile). 

ASllOl (or Sersic 159-03) is the most gas-rich system in our sample. Its enclosed gas fraction at O.lrsoo ('^0.06) is 
almost 3 times the average of other groups in the sample. Its gas fraction at r^oo (0-114j!o'o20; extrapolated) is comparable 
to those of 5 - 7 keV clusters. The Chandra exposure is short but the best-fit values of our temperature profile agree well 
with the 'K.MM-Newton results by Snowden et al. (2008). The derived C500 in this work (5.051^^ 34) is consistent with the 
result by Piffaretti et al. (2005), 4.33±0.51. 

A1991 was also studied by V05 and V06. Our derived system properties at rsoo and r25oo agree very well with those 
of V05 and V06. The temperature decline at large radii in this system was also found from the XMM-Newton data by 
Snowden et al. (2008) with consistent values. 

A2462 hosts a small corona (with a radius of ^ 4 kpc) at the center, without a large cool core. This is common for 
BCGs (Sun et al. 2007; also see §8.2). The Chandra surface brightness profile shows a significant break at ^ 180 kpc, 
which is about the size of the central radio source from the NVSS image. The central radio source in A2462 is the most 
luminous one in our sample. It may have heated the group core, as shown by the high temperature and entropy beyond 
the central corona. 

RXJ 1159+5531 has been studied by V05, V06 and G07. Our temperature profile agrees well with both V05 and 
G07. The derived gas fraction and total mass at 7^2500 are well consistent with those from V06. rsoo is also consistent. Our 
C500 (2.95^0 go) hes between the results from V06 (1.70±0.29) and from G07 (5.6±1.5). Our r^oo and M500 are consistent 
with G07's. ' 

A2550 is in a large filamentary structure that connects with A2554 (z — 0.111, 17.5' on the northeast), while A2556 
(z = 0.087) is 21' to the east. Based on the derived surface brightness profiles centered on each system, we exclude regions 
within 13.7' and 14.1' (in radius) of A2554 and A2556 respectively, which are about 1.4 times rsoo of each system. There 
is also an X-ray clump south of A2550's core that is excised. 
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Table 1 

The group sample and the Chandra observations 



Group 






ObsID" 


Date 


Exposure*^ 


PSPC 




i'1.4GHz^ 


NGC 1550 


0.0124 


51.4 


31S6 


2002-01-08 


9.7 




11.29 


21.72 








3187 


2002-01-08 


9.5 














5800 


2005-10-22 


43.0 (44.0) 














5801 


2005-10-24 


44.0 (44.3) 








NGC 3402 


0.0153 


63.5 


3243 


2002-11-05 


22.9 (29.5) 




11.40 


22.18 


A262 


0.0163 


67.8 


2215 


2001-08-03 


28.7 (28.7) 


7.6 


11.60 


22.56 








7921 


2006-11-20 


110.5 (110.5) 








NGC 383 (3C 31) 


0.0170 


70.8 


2147 


2000-11-06 


41.0 (44.3) 


24.0 


11.67 


24.46 


3C 449 


0.0171 


71.2 


4057 


2003-09-18 


14.3 (26.4) 


8.9 


11.13 


24.36 


NGC 533 


0.0185 


77.1 


2880 


2002-07-28 


29.7 (37.1) 


11.5 


11.76 


22.31 


NGC 741 


0.0185 


77.3 


2223 


2001-01-28 


28.7 (30.2) 


12.4 


11.82 


23.85 


MKW4 


0.0200 


83.4 


3234 


2002-11-24 


27.6 (29.7) 


9.4 


11.81 


22.15 


A3581 


0.0230 


96.2 


1650 


2001-06-07 


7.2 (7.2) 




11.51 


23.85 


NGC 5129 


0.0230 


96.2 


6944 


2006-04-13 


20.5 (20.5) 


5.5 


11.63 


21.90 








7325 


2006-05-14 


25.6 (25.6) 








NGC 1132 


0.0233 


97.5 


3576 


2003-11-16 


20.5 (37.1) 




11.64 


21.79 


UGC 2755 


0.0245 


102 


2189 


2001-02-07 


7.2 (15.6) 


16.4 


11.49 


24.26 


NGC 4325 


0.0257 


108 


3232 


2003-02-04 


25.6 (29.9) 


4.8 


11.31 


<21.32 


HCG 51 


0.0258 


108 


4989 


2004-02-15 


18.4 (18.9) 




11.47 


<21.32 








5304 


2004-02-16 


12.2 (12.7) 








3C 442A 


0.0263 


110 


5635 


2005-07-27 


26.6 




11.51, 11.40 


24.70 








6353 


2005-07-28 


13.6 














6359 


2005-10-07 


19.4 














6392 


2006-01-12 


31.7 








UGC 5088 


0.0274 


115 


3227 


2002-03-10 


33.8 




11.26 


21.04 


NGC 6338 


0.0274 


115 


4194 


2003-09-17 


39.9 


3.5 


11.73 


22.95 


NGC 4104 


0.0282 


118 


6939 


2006-02-16 


35.8 (35.8) 


13.6 


11.85 


21.60 


RBS 461 


0.0296 


124 


4182 


2003-03-11 


22.0 




11.41 


22.57 


ESQ 552-020 


0.0314 


132 


3206 


2002-10-14 


18.7 




11.92 


<21.72 


A1177 


0.0316 


133 


6940 


2006-12-27 


32.8 (33.5) 




11.73 


<21.50 


IC 1262 


0.0326 


138 


2018 


2001-08-23 


25.6 (30.5) 




11.49 


23.22 








6949 


2006-04-17 


38.1 














7321 


2006-04-19 


36.8 














7322 


2006-04-22 


37.4 








NGC 6269 


0.0348 


147 


4972 


2003-12-29 


38.6 


10.1 


11.94 


23.11 


ESQ 306-017 


0.0358 


151 


3188 


2002-03-08 


13.6 




11.91 


22.58 








3189 


2002-03-09 


13.8 














914* 


2000-07-26 


50.4 








NGC 5098 


0.0368 


156 


2231 


2001-08-04 


10.0 




11.44, 11.36 


23.38 








6941 


2005-11-01 


37.9 (38.6) 








A1139 


0.0398 


169 


9387 


2008-03-28 


10.0 




11.50 


23.16 


A160 


0.0447 


190 


3219 


2002-10-18 


54.5 




11.69 


24.65 


UGC 842 


0.0452 


192 


4963 


2005-02-13 


38.9 (39.2) 




11.77 


22.24 


A2717 


0.0498 


213 


6973 


2006-08-17 


46.1 


8.5 


11.82 


24.52 








6974 


2006-04-10 


18.9 








RXJ 1022-1-3830 


0.0543 


233 


6942 


2006-10-14 


40.9 (41.4) 


10.3 


11.69, - 


<21.98 


ASllOl 


0.0564 


242 


1668 


2001-08-13 


9.2 (9.4) 


11.9 


11.89 


24.27 


ESQ 351-021 


0.0571 


245 


5784 


2005-04-24 


34.8 (35.1) 




11.95 


22.71 


A3880 


0.0581 


250 


5798 


2004-12-23 


19.7 




11.93 


24.22 


A1991 


0.0587 


253 


3193 


2002-12-16 


36.8 (36.3) 


20.0 


11.86 


23.47 


A1275 


0.0637 


275 


6945 


2006-02-05 


49.1 (48.4) 




11.31 


<22.13 


A2092 


0.0669 


290 


9384 


2007-11-13 


9.7 




11.44 


<22.17 


RXJ 1206-0744 


0.0680 


295 


9388 


2007-11-15 


10.0 




11.75 


<22.19 


A1238 


0.0720 


313 


4991 


2004-03-28 


18.2 




11.70 


24.44 


A744 


0.0729 


317 


6947 


2006-10-22 


36.3 




11.80 


22.24 


A2462 


0.0733 


319 


4159 


2002-11-19 


37.9 (38.6) 


3.4 


11.82 


25.36 


RXJ 1159-1-5531 


0.0808 


354 


4964 


2004-02-11 


69.6 (74.2) 




11.97 


<22.34 


A1692 


0.0848 


372 


4990 


2004-08-12 


21.5 




11.88 


22.33 








6930* 


2006-03-06 


76.0 














7289* 


2006-03-09 


75.0 








A2550 


0.122 


550 


2225 


2001-09-03 


57.3 (58.8) 




11.62 


23.18 



"The group redshift is extracted from NASA/IPAC Extragalactic Database (NED) 
''The luminosity distance of the group derived from its redshift 

°The ObsIDs with * are observations that happened to be close to the interested groups. We used these observations to constrain the local soft 
CXB. 

''Effective exposure after excluding time intervals of background flares. For observations with ACIS-S3 at the optical axis, two exposure values 

are listed, with the one in the brackets as the exposure for the FI chips. 
'^ROSAT PSPC clean exposure time (in ksec) if the pointed PSPC observations exist. 

■'^2MASS Ks band luminosity of the cD galaxy as shown as log(LKs/L(v). Afx;;) = 3.39 mag. There are three groups (3C 442A, NGC 5098 
and RXJ 1022-1-3830) with two BCGs at the center. The Ks band magnitude of one central galaxy in RXJ 1022+3830 with two BCGs at the 
center. The Ks band magnitude of one central galaxy in RXJ 1022+3830 is unknown. 

^1.4 GHz luminosity of the cD galaxy as shown as log(Li.4GHz/W Hz~i) from the NRAO VLA Sky Survey (NVSS) or the Sydney University 
Molonglo Sky Survey (SUMSS), assuming a spectral index of -0.8, unless it can be derived from NED. 



22 



SUN ET AL. 



Table 2 

Absorption, radial range of the analysis and the local CXB 



Group 


(10^° cm-2) 


R45'' 


outermost bin^ 
(%) 


d 

det ,spe 
(kpc) 


rdet,sur 
(kpc) 


fc^i,,,,' .fo 
(ke\-) 


47-1.21kcV^ 


f 1> 

J2-8kcV 


/2 — 8kcV, expected (/limit)^ 


NGC 1550 


12.5 (10.0) 


125 


11 


355 


364 


-4-/1 l\A 

0-2olo;o7 


XI 1 u 

*-5-0.4 


1^-5-4.5' ^-4-4.0 


~12.6 (~2.1), ~8.5 (~0.56) 


NGC 3402 


11.0 (4.0) 


119 


31 


239 


239 


(0.25) 


50-0. 6 


6 5 + 2-5 
0-5-3.5 


~6.8 (~ 0.27) 


A262 


8.1 (5.8) 


152 


19 


363 


387 (800) 


0.27±0.03 


2 7+0-6 

■'■'-0.5 


in+4-0 Q d+3-4 
10_4.5. =-4-4.0 


~9.4 (0.51-1.1), ~7.7 (0.31-0.65) 


NGC 383 


5.3 


100 


21 


289 


350 (500) 


(0.25) 


„ „-t-0.4 
^■'^-0.9 


9 3 + 3-2 
■'-''-2.6 


~8.4 (0.47 - 0.61) 


3C 449 


13.3 (9.0) 


161 


32 


230 


230 (380) 


(0.26) 


6 7+0-8 
"-'-0.3 


10 0+50 


~8.4 (~0.47) 


NGC 533 


5.7 (3.0) 


103 


26 


238 


340 (380) 


0.31±0.05 


3 9+0-9 


7 1+2.5 
'-^--2.2 


~ 7.2 (~ 0.33) 


NGC 741 


5.9 (4.3) 


100 


33 


271 


360 (420) 


p 23+0-05 
O'''^-0.03 


^-^-0.3 


6 7+3-0 
"-'-3.5 


~ 7.5 (~ 0.36) 


MKW4 
A3581 


3.1 (1.8) 
5.8 (4.5) 


117 

296 


10 
25 


490 
322 


490 (720) 
450 


(0.25) 

27+0 03 
-0.02 


4 7+0.6 
*-'-0.3 
19 c:-F0.4 
^'^-5-0.2 


7 5 + 3-7 
'-"-3.9 

10.6+5-3 
— 5.5 


^8.4 (0.34 - 0.78) 
~10.1 (~0.98) 


NGC 5129 


5.5 (1.7) 


298 


45 


214 


222 (270) 


29 + 0.02 
-0.01 


12 3+0-'' 


7 0+3-0 
'-0-4.8 


~7.0 (0.17-0.41) 


NGC 1132 


7.8 (5.5) 


92 


33 


284 


310 


0.28±0.04 


^ "-0.3 


7 3+2-9 
- —2 5 


~6.5 (~0.22) 


UGC 2755 


13.6 


93 


52 


193 


193 (205) 


26 + 06 
"■'"'-0.04 


„ „+o.e 

^3-0.5 


10 4+30 
o- -3.7 


~ 8.7 (0.55-0.73) 


NGC 4325 


2.4 


143 


23 


232 


232 (250) 


(0.25) 


4 5 + 1-2 


«-li5.0 


~ 7.0 (~ 0.30) 


HCG 51 


5.3 (1.1) 


118 


30 


329 


360 


(0.25) 


2 2-1-0.6 

^-^^-0.3 


o-sia.o 


~ 7.2 (~ 0.33) 


3C 442A 


6.1 (4.8) 


85 


68 


396 


530 


(0.25) 


2 0+0.5 
^-0-0.6 


6.6±3.4 


~7.6 (0.25-0.7) 


UGC 5088 


1.2 


90 


63 


317 


320 


(0.25) 


2 5+0-3 
^-5-0.4 




~7.0 (~0.27) 


NGC 6338 


2.3 


133 


43 


349 


360 (510) 


(0.25) 


2 (1-1-0.5 
^-0-0.6 


9.2+4-1 


~7.0 (~0.31) 


NGC 4104 


4.4 (1.8) 


113 


33 


407 


420 (550) 


23+0 0^ 
"■^''-0.05 


2 7+0-4 
^-'-0.2 


«-it|:8 


~7.0 (~0.30) 


RBS 461 


17.7 (15.0) 


65 


54 


359 


515 


(0.25) 


3 6+0-3 


io-oi|:? 


~8.6 (~ 0.56) 


ESO 552-020 


3.9 


119 


51 


380 


480 


(0.25) 


2 4 + 0.4 
^■*-0.7 




~7.8 (~ 0.46) 


A1177 


4.6 (1.1) 


103 


35 


418 


420 


0.15±0.02-i 


2 n+O-"^ 


fi c+2.8 


~ 7.3 (~ 0.34) 


IC 1262 


3.4 (1.8) 


162 


74 


375 


390 


0.32±0.04 


5 4+0-4 


^■^±1:1 


6.0 - 7.3 (0.15 - 0.31) 


NGC 6269 


5.2 


106 


48 


439 


510 (740) 


(0.25) 


2.5±0.4 


»-2l3l 


~ 7.0 (~ 0.30) 


ESO 306-017 
NGC 5098 


3.0 
6.8 (1.3) 


105 
120 


77 
81 


484 
421 


490 
480 


23 + 0-04 
o-'=^-0.03 
„ 25+0 05 


4 9+0-0 

2 6+0-4 
''■°-0.5 


11 o+^-o 

11-0-2.6 

5.9+33-7, 9.5+t;| 


~ 9.7 (~ 0.85) 
~ 10.1 (~ 0.95), ~ 6.9 (~ 0.26) 


A1139 


3.2 


93 


94 


318 


370 


(0.25) 


2 6+0-4 




~ 9.6 (~ 0.87) 


A160 


4.0 


100 


62 


507 


640 


(0.25) 


3.2±0.5 


7-8ii:9 


'v^ 6.7 ('v^ 0.22) 


UGC 842 


4.0 


99 


25 


410 


510 


0.26±0.03 


4 4+0.5 
■ -0.6 


7 c+3.5 


'^6.8 ('^ 0.22) 


A2717 


1.2 


122 


37 


730 


800 (760) 


31+0 07 
"■■^ — 0.06 


3.6±0.6 




~ 6.3 (~ 0.20), ~ 8.6 (~ 0.58) 


RXJ 1022-1-3830 


3.3 (1.5) 


116 


32 


548 


610 (670) 


0.30±0.06 


2.5±0.6 




~ 7.0 (~ 0.27) 


ASllOl 


1.1 


130 


23 


568 


650 (850) 


0.28±0.07 


3.5±0.9 


^■^■"-5.0 

^■-±1:1 


~10.1 (~ 0.95) 


ESO 351-021 


2.4 


129 


34 


511 


515 


n ,„-l-0.07 


2 5 + 0-4 


~6.8 (~ 0.26) 


A3880 


1.2 


149 


37 


779 


810 


(0.25) 


o 54-1.0 


9-2ii.| 


~8.0 (~ 0.43) 


A1991 


2.5 


280 


31 


655 


820 (750) 


0.33±0.02 


14.4±0.6 


6 7+2-9 
°- '-2.5 


~7.0 (~ 0.29) 


A1275 


6.2 (2.0) 


91 


42 


495 


530 


21 + 06 


2 1+0-5 
^^-^-0.4 


7 7-F3.8 
'-3.4 


~6.3 (~ 0.16) 


A2092 


2.2 


130 


79 


555 


600 


(0.25) 


„ ,4-0.3 
■^-^-0.4 


6 3+1-9 
"-3-1.8 


~ 9.2 (~ 0.79) 


RXJ 1206-0744 


2.3 


134 


92 


451 


580 


(0.25) 


2 1+0-0 


9 2+*-3 
9-''-3.0 


~ 9.4 (~ 0.84) 


A 1238 
A744 


3.5 
3.5 


103 

104 


81 
92 


673 

433 


690 

600 


„ 2,-1-0.09 

o-'^'^-o.os 

20 + 0-30 
"■'^"-0.06 


2 4+0-4 

2 7+0-4 


8 3+4-2 
°-3-3.3 
6 4+3-3 
"-*-2.6 


~ 8.2 (~ 0.49) 
~ 6.8 (~ 0.24) 


A2462 


3.1 


158 


50 


564 


650 (800) 


0.26±0.02 


7.8±0.6 


7.6±5-0 


~6.4 (~ 0.19) 


RXJ 1159 


1.2 


134 


56 


529 


560 


0.27±0.03 


2.7±0.5 


6.0+3-1 


^5.5 (0.08-0.14) 


A1692 


1.8 


158 


80 


766 


770 


26 + 05 
O-'"'-0.04 


4.2±0.5 


7 3+3-9 


~8.0 (~ 0.39) 


A2550 


1.9 


130 


75 


597 


740 


0.28±0.03 


4.8±0.7 


7 „-F3.4 
7-»-3.1 


~ 5.9 (~ 0.13) 



"The absorption column density in our analysis. If the value from our spectral analysis (§3.3) is consistent with the Galactic value from the 
Leiden/ Argentine/Bonn (LAB) HI survey (Kalberla et al. 2005), the LAB value is used. Both values are listed if they arc significantly difi^erent, 
with the number in brackets is the LAB value. 

''The ROSAT All-Sky Survey R45 flux (Snowdcn ct al. 1997), in a unit of 10^^ cts/s/arcmin^, measured from an annulus centered on the 
source. The inner radius of the annulus is 0.4 - 0.8 dcg (depending on the source size), while the outer radius is the inner radius + 0.4 deg. In 
a few cases, we have to use partial apertures to exclude the bright sources near our targets (e.g., A1692 and A2550). 

'^The fraction of the outermost radial bin for the spectral analysis covered by the Chandra data, compajred with the full annulus. The median 
is 37%. Note this fraction is always less than one because of point sources and chip gaps. 

''The outermost radius for the Charidra spectral analysis 

''The radius where X-ray surface brightness is detected at > 2a. rout/Ti„ = 1.06-1.1. Note that our estimate of local background is conservative 
so our 2<T range is smaller than V06's 3a range for A1991 and RXJ 1159+5531. The value in brackets is for PSPC if available. 
The temperature of the hotter component of the local soft CXB 

f The 0.47 - 1.21 keV observed flux of the local soft CXB (in unit of 10~^^ ergs cm~^ deg~^). The energy band is chosen to match that of 
the RASS R45 band. 

''The 2-8 keV unabsorbed flux of the unresolved hard CXB (in unit of 10~^^ ergs cni^'^ s^'^ dcg^^) 

*The expected 2-8 keV unabsorbed flux of the unresolved hard X-ray CXB (in unit of 10^^^ ergs cm^^ s^^ deg~^), estimated based on the 
limiting flux (shown in brackets, in unit of 10"^'* ergs cm""^ s~^) and from the derived average relation in K07. 

^The 0.1 keV component has zero normalization in this case. 
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Table 3 

Derived properties of groups (I: temperature, mass and gas fraction) 



Group 



rsoo" 

(keV) 



^2500 
(keV) 



''2500 
(kpo) 



(10l3 Mr, 



/gas, 2500 



Jgas, 2500-500 



NGC 1550 


1.06±0.02* 


1 


.18±0.02 


465+;^g* 


NGC 3402 


(0.74±0.03) 





.80±0.02 


(380) 


A262 


(1.94+°11) 


2, 




(644) 


NGC 383 


(1.67+°1?) 


1, 


89"^0 10 


(593) 


3C 449 


(0.97+«-»*) 


1 


.08±0.04 


(437) 


NGC 533 


(1.06+«-08) 


1. 


2;^+0.06 


(461) 


NGC 741 


(1.27±°°|) 


1, 


3^+8:08 


(510) 


MKW4 


1.58±0.09 


1, 


7^+0-05 
04 


538+^^ 


A3581 


(i««±S:J9) 


1, 


„,+o:i2 


(593) 


NGC 5129 


(0.76±0.03) 





.83±0.02 


(384) 


NGC 1132 


(0.99±0.04) 


1 


.08 + 0.03 


(442) 


UGC 2755 


(0.76±0.05) 


0, 


OO + 0.05 
— 0.05 


(384) 


NGC 4325 


(0.89±0.03) 





.97+0.03 


(418) 


HCG 51 


(1.06±0.04, 


1. 


15+0.04 


(460) 


3C 442A 


1.34±0.n4* 


1. 


61 + 0.05 
— 06 


495+;2- 


UGC 5088 


0.81±n.03* 


n. 


„, + 0.02 
— 0.03 


364^27, 


NGC 6338 


(1.92 + S-!!?) 


2 


.14 + 0.05 


(636) 


NGC 4104 


1,41 + 009. 


1. 


64+0-0| 


535*20* 


RBS 461 


(1.93±0.10) 


2 


.17+0.06 


(637) 


ESQ 552-020 


t^- '^-0.09-' 
, ,-+0.06, 
^•^'-0.07 


1 


.97+0.09 


(598) 


A1177 


1. 


4S+0-06 
■***-0.07 


550l^9. 


IC 1262 


(1.73 + 06) 


1. 


Q,+0.06 
-0.05 


(600) 


NGC 6269 




2. 


, _ + 0.08 
I^-O.IO 


„._ + 50, 
645_52 


ESO 306-017 


„ ,-+0.12, 


2. 


54 + 0.08 
""-O.IO 


690+*** 
■'■'"-30 


NGC 5098 


0.96 + 0.04 


1. 


05+0.02 
0O-0.03 


398^:^^ 


A1139 


(2.0110J3) 


2. 


,„ + 0.35 


(650) 


A160 


1 „„ + 0.10 

i-'^**-o.io 


2. 


05 + 0.07 
-0.06 


626 + 27 


UGC 842 


1 54+0-1*. 
^■°*-0.12 


1 


.78 + 0.09 


570tf5- 


A2717 


2 43+0-13 
^■^^-0.12 


2. 


+0.08 
""-0.09 


732lf2 


RXJ 1022-1-3830 


1 94+0-20 
-'•■'*-0.14 


2. 


3„+0.17 
•^"-0.12 


631 + ^2 


ASllOl 


, J..+0.14, 

-^■"^-o.ii 


2. 


+ 0.10 
"■'-0.09 


768+|»5. 


ESO 351-021 


1 14+0-07 
-'•^*-0.04 


1, 


34+0.07 
-0.03 


437+7S 
*37_38 


A3880 


2 49 + "-" 
''•*-'-0.12 


2 


.75+0.11 


7«9if3 


A1991 


, „„ + 0.10 

-^■""-o.os 


2 


.86 + 0.07 




A1275 


1 4g+0.08, 
°-0.07 


1, 


„„+0.06 
"^-0.07 


592+|7. 


A2092 


i.67l°;13. 


2, 


14+0.17 
^*-0.21 


659l«. 


RXJ 1206-0744 


(1.91±°J?, 


2, 


+0.19 
■^*-0.21 


(624) 


A1238 




2, 


OO+0.35 
^^-0.36 


(725) 


A744 




2 


.49+0.16 


(681) 


A2462 


„ 09+0.12 


2 


.62 + 0.09 


646 + 30 


RXJ 1159 


1 84+'"'-i^* 

-^■«^_0.08 


2. 


,, + 0.09 
^^-0.10 


630t^^' 


A1692 


^■°^-0.24 


3 


+0.25 
""-0.22 


658±J* 


A2550 


1.95±0.10 


2 


.05+0.08 


617+^* 



222±6 
205^00 
288j 
252 

2iiiii 

207+ 



,+35 



7+18 
-29 



214 



+20 
-14 



259 
174 



+22 



17 



+7 
9 

215 + 11 
188+11 
212±31 
275113 

2771^5 
163^1" 
288+ J3 
274 + 12 
318+2| 

297+27 
^■"-24 
+ 20 



264 
330 
232 
337 
195 



-25 
+ 20 
-40 
+ 27 
-21 
+ 16 
-14 
+ 17 
13 

298l|^ 



286 



+ 25 
-23 

g+23 

°-29 



362 
194 



+ 28 
35 



248 
272 
325 



+63 
-59 
+30 
-24 
-+13 



3273 

273+15 
2'3_13 

356±|^ 

286+9 



„ ,„+U.32, 
■^■^°_0,37 


0.097+q 009* 





.056 + 


.003 


0.149l°;^33. 


4 93+0-5" 


(10) 






0. 


,032+0^ 


008 
.014 














0. 


,064+0- 


005 




3.48l° 


.49 


(9) 






0. 


,031 + ^' 


003 




3.09lj 




(6) 






0. 


,041 + g' 


006 




2.97I1 


'60 


(5) 






0. 


.03112; 


005 




4.58 + 3 


'90 


(5) 






0. 


,026+o' 


004 




3.05I; 


'68 
.22 


(6) 


A ar.+0.71 
4.85^Q gg 


0.086+0.009 


0. 


,047+t!' 


002 


0.13410-38 


3.93+1 


16 
.78 


(7) 






0. 


,067+S- 


006 




7.43l^ 


06 
.40 


(6) 






0. 


,035 + g- 


004 




3.43+J 


.72 
.22 


(5) 






0. 


,0391°; 


^ooi 




1 77+° 


.70 
.58 


(6) 









.031+0 


.005 














0. 


o37l"o; 


014 

.007 




5.I9I2 


.68 
.36 


(7) 









.028 + 


.003 




2 11 + 
2-11-0 


.89 
.56 


(6) 


3.90 + 0-22. 
— 0.40 


0.068 + 0-006. 
— 0.003 


0. 


028 + 0- 

— 


003 
002 


0-230 + 0-210, 
— 0. 060 








1.48 + |5'36. 







.029 + 


.003 


067+"^"21. 


4.15+; 


66 
.16 


(5) 









.053 + 


.003 




5.2712 


66 
.80 


(8) 




0-069lg:gS|* 





.036 + 


.003 


137+0.043, 


4.29I J 


.58 
.12 


(7) 









.059 + 


.006 




5.45 + 1 


.32 
.16 


(9) 









.039 + 


.005 




4-2711 


.84 
.22 


(5) 


6 28+°-*** 
o-2S_o.73 


060+°^°°^» 
o.uou_Q 007 



0. 


.037+0 

"*5l"o: 


.004 
004 
010 


n nqi+° 033, 
u.oai_o 020 


5.261* 


.68 
.61 


(6) 


S 49+1-07. 
"-*-'-2.01 


076 + °-"ll* 
"-"'6-0.010 





,044l[;; 


003 
004 


„„_+0.027, 
0,OS7_o.oi5 








io.3t2;i- 


nai +0-010. 
0-081_o oil 


0. 


"52l"o: 


003 
.004 


119+0.031, 
"■'^'^"-0.026 








"^■""-0.46 


,00+0 021 
0-108_o 012 


0. 
0. 


0481;;; 

026+o' 


006 
.004 
016 
.008 


205+°^131 
0-205_o 060 


4.2911 


.46 
.06 


(7) 


7 90 + 1-"'' 
'■""-1.10 


n 08^ + 0-009 
"■"**5_o.008 


0. 


043lg; 


003 
.004 


,,,+0.038 
0-121_o.o23 


2.73IO 


70 
.60 


(8) 


- «„ + 2.60, 
5.60_j JO* 


n nt^« + 0.012 , 
"■"56_o.014 





.032 + 


.003 


„„„+0.064, 
"■"^^-0.041 


6,06 + ^3 


31 
.19 


(6) 


12-9l?:7 


0.076 + 0.010 





.053 + 


.003 


"■"""-0.025 


2,15lS 


36 
.32 


(9) 


8.ooli;3J 


„ 07, + 0.007 
"■075_o 013 


0. 


03slo- 


004 
005 


134+0.074 
u.ia4_o.o46 


4.031J 


24 
.80 


(6) 


14 1 + 5.0, 


114+0.021, 
"■-'-'*-0.020 





.085 + 


.008 


129 + 0-OSS, 
u.i2J_o 041 


5 05 + 2-37 
5-"''-1.34 


(10) 


„ ,,+ 1.80 
3 -0.90 


074 + "-"13 
"■"'*-0.018 


0. 


,036+2^ 


007 
.005 


,,Q+0.066 
"■^^"-0.050 


2.42IJ 


20 
.76 


(6) 


-■^ir.°5 


088 + 0-016 
0. 088^0 021 


0. 


,0741°: 


005 
.007 


090+°^°*'' 
u.u»u_o 022 


4.03+ J 


54 
.26 


(8) 


13.4+2;| 


094+°^"l" 
u.uy4_o 012 





.066 + 


.004 


,,-+0.038 
0-11O-0.022 


*-69l°o 


.76 
.70 


(9) 


6 90+3-00 * 
"'■""-1.67 


094+0 °l** 
u.ua4_o 031 





.069 + 


.004 


ioo+° 032, 
u.iuu_o 035 


4.I0II 


.69 
.27 


(7) 


8.95ii;«i* 


0.078+0.013* 


0. 
0. 



,0501°; 

,0531°; 
.035 + 


006 
.007 
009 

.007 
.008 


096+°^°^''* 
U.U»D_o 020 












0. 


.046lg; 


005 
.006 




3.33+ J 


66 
.04 


(6) 


8.80ti;f, 


099+"^"ll 
"-"""-0.009 





.052+0 


.003 


igo+°-058 
u.nsu_o 034 


3.41IJ 


.56 
.84 


(6) 


„ ,0 + 3-10, 
8.30_j j2 


n n«t^ + 0.007, 
"■"6S_o.012 





.042+0 


.003 


o=n+°-019» 
°-080_o 023 


2.95IJ 


.16 
.90 


(6) 


9 70+2-99 
"■70-1.91 


090 + "-"l* 
"■"""-0.020 





.043 + 


.006 


216 + "-155 
"■216-0. 087 


5.46+1 


.77 
.19 


(7) 


7 90+2-90 
'■""-1.00 


093 + " "ll 
U.UJJ_o 016 


0. 


06712; 


003 
004 


119 + 0-027 
u.i.i.j_o.038 


4 57+0 


.92 
.76 


(7) 



"The T500 in parenthesis are from the empirical relation between T1500 and T500 (§4). The rsoo in parenthesis are estimated from the M500 — T500 
relation in this work. The values with an asterisk are for tier 2 groups from extrapolation. 

''The value in parenthesis is the number of radial points between 40 kpc and rjet spe (including 40 kpc and rjet ^pe) 3-re used in the NFW fit 
(§7.3). 
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Table 4 

Derived properties of groups (II: entropy) 



Group 




-'^ looo" 


-^1500'* 


-K2500 


-'^O.lSrSOO 


-^30kpc 


NGC 1550 


(297i:») 


253±23 


228±15 


1981" 


94±3 


49±1 


NOG 3402 






(2781^1**) 


226+ if 


66±4 


34±2 


A262 






3831^^ 


3541S 


I61I? 


69±2 


NGC 383 






(495lfe) 


4481^^ 


327+ 


2141?? 


3C 449 






(217lf,) 


199+27 


i58i;i 


96±6 


NGC 533 




(3361??) 


317l« 


289l^J 


17111? 


103±6 


NGC 741 




(460t«tj 


42811^ 


388111 


243l|^ 


14611? 


MKW4 


574ire 


454tr, 


3881- 


332+^^ 


161±6 


8OI? 


A3581 




(384l«4) 


344I59 


298lt? 


124±7 


38±1 


NGC 5129 




(2621^?) 


229l| 


I93I1? 


100±6 


751? 


NGC 1132 






336ll« 


238±26 


-6 


68+? 


UGC 2755 






(260+^^) 




120111 


113l?t 


NGC 4325 






(195l?«) 


^ ' ^-36 


841?, 


39±2 


HCG 51 




(37lllfe) 


348lf„? 


303+S 


126±5 


82±4 


3C 442A 


(455l«) 


438ir3 


4311- 


41911? 


210l?^ 


72±8 


UGC 5088 




326+tl 


288l|^ 


230±23 


84+? 


5111 


NGC 6338 




(40llt«) 


071 +36 


3341^1 


239+11 


108±5 


NGC 4104 




4711- 


42III? 


36O+2I 


221II? 


14811? 


RBS 461 






(446lfe) 


370lt? 


165±8 


93±5 


ESO 552-020 




(580ll°I) 


550ir. 


460l^? 


190±9 


97l? 


A1177 






397+f, 


317+^^ 


179+" 


13811? 


IC 1262 






(489l«|) 


411+™ 


134±4 


47±2 


NGC 6269 






434ir„ 


40615^ 


272II? 


1271^ 


ESO 306-017 






5961^;^ 


475+" 

^''-'-38 


191±12 


86±6 


NGC 5098 


217±43 




205±21 


194+;? 


84+^ 


45±2 


A1139 




(486l-§) 


43^+112 
^■^-^-100 


36811? 


26111? 


A160 


477-1-86 


429tTe 


407lf« 


379+?? 


26811^ 


172l?i 


UGC 842 


(6521^^1) 


522+3° 


4591™ 


4061^3^ 


2381?^ 


I22I11 


A2717 


1022+J^^ 




673+^1 


50211? 


186±11 


75±5 


RXJ 1022+3830 


655t^,^ 


587l^,« 


5481^5? 


4881S 


2231?? 


77l? 


ASllOl 


(476lf,D 


465+;,f 


441+73 


386l|? 


I20II? 


29±1 


ESO 351-021 


309+^f 


o-|r-|-S5 


310l|^ 


2891^1 




57±4 


A3880 


624l^fe 


68ii;,«/ 


560+^^ 


420+*^ 


176±20 


50±3 


A1991 


752+^^^ 


6531?=, 


568lf, 


4631- 


I66+I3 


39±1 


A1275 


(464tf/) 


324l«« 


284I30 


231+?? 


I30III 


59±3 


A2092 


(42311!^) 


449+00 


4391^9 


401+?? 


202+?l 


200l?l 


RXJ 1206-0744 




4321-, 




354+?? 


201III 


102+?? 


A1238 




8051?- 


676i;«^ 


584111/ 


464llf 


34^+200 


A 744 




(695l=,fo) 




470IS 


21511? 


11511? 


A2462 


627t\1l 


5371™ 


488lfo 


42211? 


24711? 


1721?? 


RXJ 1159 




580l^« 


491+*3 


4101?? 


230II? 


92±6 


A 1602 


j^29 ~f" 202 




r--jr + 174 


rv(- + 126 

57o 


252 + 29 


14!) + ?? 


A2"i:.() 


1 ii 


:.:!() J"? 


l:u ' - 


?? 


111'?^ 


l(j 1 :i 



"The values in parentheses are from extrapolation (see §4). 
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Table 5 

K - T relations" 



Relation 


Ki 


a 




(kcV cm^) 




K500 - T500 (tier 1) 


286+60 


1.12+0.28 


K500 - T500 (tier 1+2) 


356+45 


0.825+0.200 


Kwoo - T50Q (tier 1+2+3) 


310+19 


0.832+0.100 


i^i500 - T500 (all groups) 


288+14 


0.790+0.078 


- ^2500 (all groups) 


230+10 


0.760+0.061 


-K'o.i5r50o - ^2500 (all groups) 


114+9 


0.778+0.124 


^30kpc - ^2500 (all groups) 


68+9 


0.459+0.229 


K500 - Tsoo (tier 1 + clusters) 


286+23 


1.08+0.05 


- T500 (tier 1+2 + clusters) 


329+25 


0.994+0.054 


-fi^iooo - T500 (tier 1+2+3 + clusters) 


303+13 


0.887+0.034 


-f^2500 - 7500 (all groups + clusters) 


252+8 


0.740+0.027 


-f^o.i5r500 - 7500 (all groups + clusters) 


137+9 


0.494+0.047 



°-E(z)^/'^-ft' = Ki (T / 1 kcV)™, where Ki is the corresponding entropy at 1 keV. The last four relations include 14 clusters from V06 and V08. 
There is little correction between -ftTsokpc ^^nd T250O) but we still list it for the studies of intrinsic scatter (Fig. 11). We used the BCES (Y|X) 
regression (Akritas & Bershady 1996) as the temperature errors are smaller than the entropy errors. 



Table 6 
M500 - Tboo relation'' 



Sample 


M3 


a 


rs 




(lO"/i-iM0) 




(/i-i Mpc) 


Tier 1 


1.17+0.21 


1.64+0.21 


0.587+0.035 


Tier 1 + clusters 


1.21+0.08 


1.68+0.04 


0.593+0.014 


Tier 1+2 


1.27+0.12 


1.67+0.15 


0.602+0.020 


Tier 1+2 + clusters 


1.26+0.07 


1.65+0.04 


0.600+0.011 


Tier 1+2 (> 1 keV) + clusters 


1.27+0.06 


1.60+0.03 


0.602+0.010 



'^E{z)Msoo = (T500 / 3 keV)", where T500 is spectroscopic temperature, is the corresponding scale in the rsoo — Tsoo relation {E{z)r5oo = 
(Tsoo / 3 keV)"/'^). The cluster sample includes fourteen T500 > 3.7 keV systems from V08. We used the BCES orthogonal regression 
(Akritas & Bershady 1996). 



Table 7 
Msoo - Vx,6oo relation'' 



Sample 


Mc 


a 




(lO^ft-iMo) 




Tier 1 


1.12+0.16 


0.605+0.050 


Tier 1 + clusters 


1.10+0.07 


0.588+0.012 


Tier 1+2 


1.12+0.10 


0.564+0.031 


Tier 1+2 + clusters 


1.14+0.05 


0.571+0.010 



°-E{z)'^/^Mzoo = Mc (>x,500 / 4x10" keV Mq)". The cluster sample includes fourteen Tsoo > 3.7 keV systems from V08. We used the BCES 
orthogonal regression (Akritas k, Bershady 1996). 
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ovn/o VIII 

from the soft CXB 



Fe L hump from 
' the group emission 



0.5 




1 2 
channel energy (keV) 



Fig. 1. — The spectra of the NGC 1550 in the outermost radial bin that is to the edge of the Chandra FOV (one from the SI chip of the 
ObsID 5800 in the solid line and another from the S2 chip of the ObsID 3187 in the dashed line). Besides the still significant iron L hump 
from the group emission, the O hump from the soft CXB is also strong in the SI spectrum, which allows a robust separation of these two 
components. 




100 200 300 400 

RASS R45 (lO"*" cnts s^^arcmin"^) 



Fig. 2. — The 0.47 - 1.21 kcV observed flux surface density of the soft CXB from the Chandra data vs. the RASS R45 flux (in RASS 
channels of 13 - 18) measured in an annulus around the target (inner radius of 0.4-0.8 deg and the outer radius is the inner radius + 0.4 deg). 
The three lines arc the expected conversions between two fluxes, with the assumed two thermal components for the soft CXB (all with Teool 
= 0.1 kcV). The solid line is for: Thot = 0.2-5 keV and NORMhot / NORMcool = 0.5. The dotted line is for the T^ot = 0.2 kcV and NORMhot 
/ NORMcool = 0.5. The dashed line is for: Thot = 0.3 keV and NORMhot / NORMcool = 1- The assumed absorption is 4x10^" cm-^. 
The total hard CXB (1.74x10^-'-' ergs s^^ cm^^ dcg^^ in the 2-8 keV band from K07) is added in the conversion. The good agreement on 
average can be seen. One should be aware that the RASS R45 flux is extracted in a much larger area surrounding the interested group range 
and some of the soft CXB may come from the time-variable SWCX emission. We also include 3C 296 (T^ot = 0.38 ± 0.03) in the high CXB 
flux end (star). It is an 1 keV group but was left out as the high local soft CXB (on NFS) prohibits deriving gas properties to r25oo- 
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Fig. 3. — The temperature profiles of 43 groups in our sample. The red lines are the reconstructed 3D temperature profile with the 1 a 
errors, derived from 1000 simulations. The black line is the best-fit projected temperature profile. The effective radius of each bin for the 
projected temperature profile is derived by weighting the projected emissivity profile with the actual spatial coverage of the each bin. r25oo 
and 7*500 are marked. The rsoo in parentheses is estimated from the M500 — T500 relation. 
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Fig. 4. — Continue from Fig. 3. 
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Fig. 5. — Continue from Fig. 4. 
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1 10 100 

Radius (kpc) 

Fig. 6. — The best constrained 3D abundanee profile in tliis sample with the best parametric fit (§3.3) and 1 cr errors from 1000 simulations. 
The 1000 simulated abundance profiles are all used to derive 1000 simulated temperature and density profiles. As errors of other abundance 
profiles in this sample are larger and there are fewer bins, our simple 3D abundance model (§3.3) always fits well beyond the central 10 kpc. 
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Fig. 7. — One example of the surface brightness profiles {Chandra + PSPC) with the best fits derived from the best-fit 3D temperature 
and abundance profiles. The Chandra profile is the upper one, while the PSPC profile is the lower one. We generated response files for each 
Chandra radial bin. Note that the density errors are derived from 1000 Monte Carlo simulations with 1000 simulated 3D temperature and 
abundance profiles folded in. The dashed and dotted lines mark r2600 and rsoo (see Fig. 4 for NGC 6269). 
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Fig. 8. — Temperature profiles in the logarithmic and linear scales of rsoo- Temperatures are scaled by T25oo- Despite the large scatter at 
small radii, the temperature profiles outside of 0.2 rsoo are generally similar. The thick solid line in the linear plot is the universal temperature 
profile (also projected) derived from the simulations in Loken et al. (2002). We simply used r25oo to replace Tq in Loken et al. (2002). Good 
agreement can be seen even though the normalization is not adjusted. The thick dashed line is a simple linear fit to the data (see §5). 




500 

Fig. 9. — Mean temperature profile of groups (black circles) and the 1-it scatter in dotted lines. The solid line is the best-fit from equ. 6. 
The dashed line is the mean temperature profile of 1-3 keV systems from Borgani et al. (2004) simulations. The data points in red triangles 
are the mean temperature profile from LM08 on 48 kT > 3.3 keV clusters at 2 = 0.1 - 0.3. Note the mean temperature Tm defined in LM08 
is computed by fitting the profile with a constant after excluding the central 0.1 rigo region. It should be smaller than T2500 as T2500 is 
emission- weighted within ^2500- Nevertheless, it is clear that the group temperature profiles are more peaky than those of clusters around the 
center. 
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Fig. 10. — The entropy profiles scaled by Xsoo.adi (§6). The dashed line represents the basehne entropy profile derived by VKB05 with a 
power index of 1.1. The observed entropy profiles all lie above the baseline. 
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Fig. 11. — Entropy values at rsoo, J'looo, risoo, ^"2500, 0.15 rsoo and 30 kpc vs. system temperature (T500 or Tgsoo). The open data points 
(in the upper three panels) are based on extrapolation of the temperature and density profiles (see §4 for details). The solid lines are the 
best-fits to our data from the BCES (Y|X) estimator (Table 5), while the dashed lines axe the best-fits from the NKV07 simulations. The 
entropy excess above the NKV07 simulations (with cooling and SF) is significant at r2600, while the agreement is better at larger radii. The 
dotted line in the A'soo plot represents the base-line entropy by VKB05 (or 1.40 J^soo adi> §6). Note that the baseline has a slope of 1.1 rather 
than 1.0 as the M500 — T500 relation used (the fourth row of Table 6) has a slope of 1.65 (rather than 1.5). It is also clear that the entropy 
values at 0.15 rsoo and 30 kpc show large intrinsic scatter. The measured intrinsic scatter in the K — T relations decreases with radius and 
stays the same from rgsoo to rsoo at ~ 10%. 
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kTsoo (keV) kJ^,, (keV) 

Fig. 12. — K — T500 relations at rsoo, ''lOOOi ''2500 and 0.15 rsoo for groups in our sample and 14 clusters from V08. The solid lines are 
the best-fits of the data from the BCES (Y|X) estimator (Table 5), while the dashed lines are the best-fits from the NKV07 simulations. The 
dotted lines in the K^qq and -^^2500 plots represent the base-line entropy by VKB05. The best-fit slopes from observations are also shown. The 
agreement between observations and the NKV07 simulations becomes better with increasing radius. At r2r)()(). the observed entropy values 
are on average 56% - 22% higher than those from the NKV07 simulations at 0.8 - 2.5 keV. At riooo, the difference is 18% - 8% from 0.8 - 
10 keV. At rsoo, the NKV07 line is basically the same as our best-fit, which also has a slope expected from the self-similar relation (1.0). 
The K().i5r500 — Tsoo relation has significant scatter. Almost all clusters in the V08 sample have dense cool cores. Inclusion of non-cool-core 
clusters may steepen the K — T relation at 0.15 rsoo. 
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Fig. 13. — The ratios of the observed entropy values to the entropy baseline from VKB05 at rsoo, riooo and r2500. Fourteen clusters from 
V06 and V08 arc included for comparison. The open data points are from extrapolation. We also show weighted means at each radius for 
the group and cluster samples. The observed entropy values are always larger than or comparable to the baseline at all radii, but the average 
ratios decrease with radius for both clusters and groups. The decrease is more dramatic in groups. 
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Fig. 14.— The entropy slopes between risoo - rdet,spo, f2500 - ?'1500 (or r-det,spe), 0.15 rsoo - r2500 and 30 kpc - 0.15 rsoo vs. Tsoo- Beyond 
0.15 rsoo, the slopes are always shallower than 1.1 and the weighted averages are all around 0.7. 
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Fig. 15. — Afsoo — T500 (left panel) and M500 — ^x,500 (right panel) relations, combining the results from this work and V08. The open data 
points are tier 2 groups. The solid lines are the relations for the real mass in the NKV07 simulations, while the dotted lines are the relations 
for the mass derived under the assumption of hydrostatic equilibrium in the NKV07 simulations. The dashed lines are the best-fit relations 
from V08. The M500 — Tsoo relation can be well described by a power law down to at least M500 of 2 xlO^'^ h~^ Mq, although the HSE 
mass may be systematically lower than the real mass. The M500 — 500 relation has a smaller scatter and the agreement with the NKV07 

simulations is much better. Note as Ix.soo °c the /i-dependence of M500 should be h^/^ for the self-similar relation, M500 oc ^x^^oO' 

We still use the dependence to directly compare with the left panel 
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Fig. 16. — The density gradient at rsoo in terms of /3efF (see the appendix of V06) for tier 1 and 2 groups that M500 is derived. Most groups 
have /3eff of 0.55 - 0.75. 
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Fig. 17. — The enclosed gas fraction within r25oo vs. T500 (groups + 14 clusters from V06 and V08). Tlie solid line is the BCES fit to the 
group sample ro.0355±0.0018 (Tsoo/l keV)" while the dotted line is the BCES fit to the group + cluster sample (0.0347±0.0016 

(T500/I keV)°-^09±0 034). The intrinsic scatter on the /gas ,2500 — T500 relation is 22%. Two dashed lines enclose the 1 a region of the universal 
baryon fraction derived from the WMAP 5-year data combined with the data of the Type la supernovae and the Baryon Acoustic Oscillations 
(0.1669±0.0063, Komatsu et al. 2008). 
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Fig. 18. — The enclosed gas fraction within r25()o is correlated with the entropy at 0.15 rsoo, after the temperature dependence on both 
variables are removed (right panel). The slope is about -0.7. The intrinsic scatter on the scaled /gas, 2500 and the scaled ii'o.isrsoo is 11% 
and 14% respectively, compared with 22% and 29% intrinsic scatter in their relations to temperature. This is primarily driven by density - 
density correlation, but quantitatively shows their connection. 
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Fig. 19. — The enclosed gas fraction within rsoo (left panel) and between r25oo - ^'500 (right panel) vs. T^oo. The open data points are 
tier 2 groups. The solid lines are fits to all 23 groups and 14 V08 clusters, 0.0708±0.0046 (Tsoo/l keV)"-^^"*" "'"' , while the dashed line in 
the left panel is the fit excluding 12 tier 2 groups, 0.0776±0.0057 (T500/I keV)"-i58±0.047^ ^he listed power law index can be compared with 
that in NKV07 simulations (0.152). Two horizontal dashed lines enclose the 1 a region of the universal baryon fraction (see the caption of 
Fig. 17). We also estimate the total baryon fraction within rsoo by adding the /gas. 500 ~ ^500 relation from 23 groups and 14 V08 clusters 
and the relation for stellar fraction from Lin et al. (2003) and Gonzalez et al. (2007), as shown by the two dotted lines. The upper one uses 
the Gonzalez et al. relation, while the lower one uses the Lin et al. relation (see §7.4, note that the Gonzalez et al. relation includes the 
intracluster light). It appears that gas fraction between r25oo - ''500 has no temperature dependence (/gas ~ 0.12) on average, although there 
is still scatter. The solid line is the fit to all 23 groups and 14 clusters with a slope of 0.018±0.068. The dashed-dotted line is the expected 
average /gaB,2500-B00 expected from the best-fit /gas, 500 — T500 and /gas, 2500 — Tsoo scaling relations in this work (see §7.2). 
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Fig. 21. — C500 — Afsoo relation from this work. The filled and open circles are for tier 1 and 2 groups respectively. The open triangles are 
for tier 3 and 4 groups with mass estimated from the M^oo — Thoo relation derived in this work. The solid line is the median relation from 
the model of Bullock et al. (2001) with parameters F=0.001 and K=2.8, while the two dotted lines enclose the 1<t region. We use: <T8=0.817, 
f^M=0-279, ns=0.96 (tilt) from Komatsu et al. (2008). Note that affects the predicted C500 significantly (see Buote et al. 2007). The 
dashed line is the best-fit from GOT. Our results show no significant mass dependence of C500 in this narrow mass range and are generally 
consistent with the prediction under the current value of as- The weighted mean of C500 is 4.2. 
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Fig. 22. — 2MASS Ks band luminosity of the cD galaxy vs. A/soo- The filled and open circles are for tier 1 and 2 groups respectively. 
The open triangles are for tier 3 and 4 groups with mass estimated from the M500 — T^oo relation derived in this work. The red points 
are fossil groups identified in this work. There is a general trend that more massive groups host more massive central galaxies. The dotted 
line represents a constant A/500 /iKs,cD line. For the same system mass, the fossil groups host more luminous (or more massive) cDs than 
non- fossil groups. 




Fig. 23. — Two groups with a significant entropy bump: UGC 2755 and 3C 449. Both hosts an active FR I radio source (see §8.1). The 
dotted lines enclose the Icr error region, while the dashed lines represent the baseline entropy profile (estimated from the group's mass) with 
a slope of 1.1. 
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Fig. 24. — The spectral residual after the period D stowed background subtracted by the period E stowed background (in black). Their 9.5 
- 12 keV fluxes are scaled to be the same. The residual is flat and is only ~ 6%. The red data points are the residual (~ 0%) if the period D 
stowed background is scaled up by 5.7%. Clearly after this special scaling, the 9.5 - 12 keV fluxes of two background are not the same, which 
implies a small spectral change for the BI PB from the period D to E. 



